Hypothesis of two-dimensional stripe arrangement and its implications for the 
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The hypothesis that holes doped into high- Tc cuprate superconductors organize themselves in 
two-dimensional (2D) array of diagonal stripes is discussed, and, on the basis of this hypothesis, 
a new microscopic model of superconductivity is proposed and solved. The model describes two 
kinds of hole states localized either inside the stripes or in the antiferromagnetic domains between 
the stripes. The characteristic energy difference between these two kinds of states is identified with 
the pseudogap. The onset of superconductivity is caused by the interaction, which is assumed to 
be mediated by the transverse fluctuations of stripes. The superconducting (SC) order parameter 
predicted by the model has two components, whose quantum phases exhibit a complex dependence 
on the the center-of-mass coordinate. The model predictions for the tunneling characteristics and 
for the dependence of the critical temperature (Tc) on the superfiuid density show good quantitative 
agreement with a number of experiments. The model, in particular, predicts that the SC peaks in 
the tunneling spectra are asymmetric, only when A/Tc > 4, where A is the SC gap. It is also 
proposed that, at least in some high-Tc cuprates, there exist two different superconducting states 
corresponding to the same doping concentration and the same critical temperature. Finally, the 
checkerboard pattern in the local density of states observed by scanning tunneling microscopy in 
Bi2Sr2CaCu2 0g4.5 is interpreted as coming from the states localized around the centers of stripe 
elements forming the 2D superstructure. 

PACS numbers: 74.72.-h, 74.20.Mn 
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I. INTRODUCTION 

At present, the school of thought, which stipulates 
the importance of local inhomogeneous structures called 
"stripes" for the physics of high- Tc cuprate superconduc- 
tors, is dominated by the view that stripes form a one- 
dimensional(lD) array, where each of them runs parallel 
to one of the principal lattice directionsiiSi^iiiSiSiLSiSiiSiii. 
This picture, however, entails a number of difficulties as- 
sociated with the stripe geometry. In particular, in the 
presence of stripes parallel to the principal lattice direc- 
tions, it is difficult to explain why the "nodal" quasi- 
particles having momentum directed along lattice diago- 
nals (i.e. at 45 degrees with respect to stripes), remain 
least gapped with the onset of both the pseudogap phe- 
nomenon and superconductivity. 

In this work, I propose a microscopic model, which is 
based on the hypothesis that, at those doping concen- 
trations, where superconductivity is observed in high- Tc 
cuprates, the holes organize themselves into a two- 
dimensional (2D) array of diagonal stripes. The 2D stripe 
superstructure does not incur the geometrical difficulties 
associated with ID stripe arrays. This superstructure 
has been mentioned in the literatureii^iSiiiiiSiH (some- 
times under the name of "grid" or "checkerboard" ) , but a 
number of experimental and theoretical arguments have 
been put forward against its existence. However, as I 
discuss later (in Section [fflj, the 2D stripe scenario has 
not been yet ruled out entirely. At the same time, this 
scenario has never been analyzed persistently enough, in 
part, because no theoretical model has been put forward, 
which would relate the 2D stripe superstructure to super- 



conductivity. 

The model proposed in this work reconciles the 2D 
stripe geometry with superconductivity by demonstrat- 
ing, that the superconductivity can be carried by states 
localized in the 2D stripe background. The interaction, 
which, in the model, leads to the superconducting (SC) 
transition, is, presumably, mediated by the transverse 
fluctuations of stripes. 

The picture emerging in the framework of the 2D di- 
agonal stripe hypothesis also offers a very simple inter- 
pretation of the pseudogap phenomenorki^ii^iiSiii, includ- 
ing its role in the superconductivity of cuprates. Other 
experimental facts to be interpreted in this work are: 
(i) Quasiparticle coherence in fc-space, which emerges 
only below the SC transition; (ii) The asymmetry in 
the tunneling density of states; (iii) Linear density of 
states in the vicinity of the chemical potential; (iii) The 
checkerboard pattern in the local density of states in 
Bi2Sr2CaCu208+5 (Bi-2212) observed by scanning tun- 
neling microscopy (STM)i24i2iSS42i (iv) Low superfiuid 
density and the universal dependence thereof on the crit- 
ical temperature. 

The SC order parameter obtained in this work has 
complex two-component structure, which can not be de- 
scribed as either s-wave or d-wave or the combination of 
two. The distinctive unconventional feature of this order 
parameter is the non-trivial symmetry with respect to 
spatial translations, which includes the sign change of at 
least one of the two components. 

Reviewing the relevant literature, it should be noted 
that the general idea of superconductivity carried by 
localized states has been discussed in the past in the 
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context of various physical systems including high- 
Tc cuprates22i2i2iiS^. At a different level, this 
work also contains parallels with several theoretical 
proposalaSLSSiSS, which involve two-component scenarios 
for high-Tc cuprates. Since the 2D stripe superstructure 
can be viewed as a collection of nanodomains having 
different electronic properties, the present work can be 
linked to a more general class of ideas stipulating some 
kind of phase separation in cuprates (see, e.g., Ref.'^"). 
If only the charge ordering associated with the 2D stripe 
superstructure is considered, then it exhibits certain sim- 
ilarities with Wigner crystals advocated in some of more 
recent proposals^iiS^. 

The rest of the paper is organized as follows: In Sec- 
tion m the hypothesis of 2D diagonal stripe configura- 
tion is formulated. In Section IIIII the arguments in fa- 
vor and against the existence of such a configuration are 
discussed. Possible dynamic properties of this configu- 
ration are briefly analysed in Section Hvl In Section Ivl 
the properties of hole excitations in the presence of the 
2D stripe configuration are described, and the pseudo- 
gap phenomenon is identified. In Section IVII a model 
describing the hole excitations in the stripe background 
is formulated and solved in the mean-field approximation. 
In Section lVlII the realistic features of the model and the 
resulting phase diagram are discussed. The model pre- 
dictions are compared with experiments in Section [VIIII 

This paper is quite long, in part, because some of the 
model predictions tested in Section [VIIII require detailed 
calculations. In the first reading, one can, therefore, re- 
view all the figures in the theoretical sections, read Sec- 
tion and then proceed with reading Section IVlIII 



II. TWO DIMENSIONAL CONFIGURATION OF 
DIAGONAL STRIPES 

In this section, I introduce the basic assumption of the 
present work, namely, the two-dimensional configuration 
of energetically deep and spatially narrow stripes. 

I assume that, at sufficiently high doping concentra- 
tions, high-Tc cuprates find it energetically favorable to 
organize the spins of copper atoms into the background 
of antiferromagnetic (AF) domains as shown in Fig.^a). 
Such a background creates an effective potential for the 
holes with minimum at the boundaries between those 
domains^. Indeed, if one considers only the nearest 
neighbor exchange interaction between spins, then plac- 
ing a hole in the middle of an AF domain would cost 
energy 4J, whereas at the domain boundary it costs J. 
(Here J is the nearest neighbor exchange coupling con- 
stant.) The hole kinetic energy at the boundary is also 
lower than inside an AF domain, because in the latter 
case, the hole cannot hop to the neighboring sites with- 
out increasing the exchange energy of the system, while, 
in the former case, it can. I further assume that the gain 
in the exchange and the kinetic energies outweighs the 
loss in the Coulomb energy (associated with the repulsion 



between holes). Therefore, holes fill the domain bound- 
aries and thus form stripes (see Figs.^b) and^c)). It 
is not important for the present work whether stripes are 
centered on copper atoms, as shown in Figs. ^b,c), or 
on oxygen atoms (i.e. on the "bonds" between copper 
sites). 

The AF domains formed between the stripes fall in two 
groups, which can be distinguished by the value of an AF 
index 77 (-1-1 or —1) representing the sign of the AF order 
parameter within a given domain. It is easy to see in 
Figs, n^a-c), that the AF indices 77 of two neighboring 
AF domains always have opposite signs. 

A few mathematical facts about this kind of superlat- 
tice: 

The spin periodicity of such a structure along each of 
the two diagonal directions is 2/, where I is the side length 
of a single AF domain. 

For such a structure, the main splitting of the mag- 
netic inelastic neutron scattering (INS) peak around the 
AF wave- vector ( — , — ) is expected to be four-fold as 
shown in Fig. |2 (Here ao is the lattice period.) The 
characteristic wave vector qp of this splitting is frequently 
parametrized as 



gp = — 0, 



where, in the present case, 
5 = 



(1) 



(2) 



Theoretically, there should also exist other peaks corre- 
sponding to the higher order Fourier harmonics of the 
stripe superstructure. 

Neglecting the intersections between stripes, the frac- 
tion of lattice sites lying on the stripes is 
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The assumption of energetically deep stripes made earlier 
implies that most holes are located inside the stripes. In 
this case, the fraction of stripe sites occupied by holes is 
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where Xd is the dimensionless doping concentration. 

Yamada et alj^ have discovered experimentally that 
in underdoped La2-a;Sra;Cu04 (LSCO) 
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which means that c « ^ , or the hole content of diagonal 
stripes can be set empirically at 1 hole per two stripe 
sites (the same as 1 hole per length 2aoV2). 

Below, I will use term "stripe element" to refer to the 
part of a stripe, which constitutes a side of one AF do- 
main. On the basis of empirical relation (O, I can esti- 
mate the number of lattice sites in one stripe element as: 
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FIG. 1: (color online) 2D configuration of diagonal stripes: 
(a) Stripes in spin structure isolating two AF domains with 
opposite sign of AF vector, (b) The same as (a) but with 
holes (filled circles) occupying every second site along the spin 
stripes. Spins on the remaining stripe sites are circled to guide 
the eye. (c) The same as (b) but with more stripe supercells 
shown. 
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FIG. 2: Four-fold splitting of magnetic inelastic neutron scat- 
tering peak. 
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Doping concentration Xd — 1/8, taken as an example, 
then corresponds to I = 4aoV2- 



III. ARGUMENTS IN FAVOR AND AGAINST 
THE EXISTENCE OF THE 2D DIAGONAL 
STRIPE STRUCTURE 

If stripes exist, then the primary evidence for their 
geometrical properties comes from the four-fold split- 
ting of the magnetic INS peak (usually called the (tt, tt) 
peak), which was described in the previous Section and 
shown in Fig. |21 This splitting was observed in the 
underdoped compounds of LSCO^ and YBa2Cu306+x 
(YBC0)^S*^Si2i families of high-Tc superconductors. 

A straightforward interpretation of this peak pat- 
tern is that the antiferromagnetic spin structure 
cos(^a::) cos(^y) is modulated along both diagonal di- 
rections by function cos[^(a:; -I- y)] cos[^(a; — y)]. (Co- 
ordinates X and y correspond, respectively, to the hor- 
izontal and vertical axes in Fig.l(a-c).) The mod- 
ulation function can then be expanded as the sum 

— hence the 



four-fold splitting of the main peak. 

However, since the early indications of the stripe na- 
ture of the (tt, tt) peak splitting^, most of the "stripe com- 
munity" has opted for an interpretation of experiments 
in terms of the superposition of two two-peak splittings. 
The two kinds of splitting were described as coming from 
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two kinds of spatially separated domains — each repre- 
senting a one-dimensional array of stripes running along 
one of the principal lattice directions. 

There, indeed, exist a number of theoriesi^SiSi^iSiSii sug- 
gesting that stripes extending along the principal lattice 
directions are more favorable energetically than diago- 
nal stripes. (The 2D diagonal stripe superstructure has 
been explicitly considered, e.g., in Refsi^*iS.) This is, 
however, a delicate energetic balance, which should be 
sensitive to numerous factors, not all of which are taken 
into account by the above mentioned theories. For ex- 
ample, the interaction with the crystal lattice and the 
long-range Coulomb interaction are, typically, neglected 
in the numerical studies, even though the energy asso- 
ciated with each of these two interactions can change 
the outcome of the competition between different stripe 
configurations. Since the chances of bringing the stripe 
energetics under the full control of first principles calcu- 
lations are quite slim, the choice between different stripe 
configurations (including the absence of stripes) should 
be made, eventually, on the basis of experiments. 

On the experimental side, one can find three main ob- 
servations suggesting the one-dimensional nature of the 
stripe pattern^ ». 

First, in other materials structurally similar to high- Tc 
superconducting cuprates, such as, e.g., insulating 
Lai.gsSro.osCuO^^- and nickelate La.2-xSrxNi04+jj^^i^, 
there exists direct experimental evidence of the ID na- 
ture of stripe modulations. 

Second, the two-dimensional spin structure modula- 
tion of the form cos[^(a; + y)] cos[^(x — y)] should in- 
duce an effective potential for electric charge roughly of 



the form ^cos[^{x + y)] cos[^{x - y)]j . If the spin 

modulation is weak, then it follows from the Landau the- 
ory of the second order phase transitiong^^'^^, that the 
peak structure for the charge inelastic neutron scatter- 
ing corresponding to the above potential should exhibit 
four main peaks around (0,0). The orientation of these 
peaks should be rotated by 45 degrees with respect to 
the splitting of the magnetic peak, and characterized by 
the separation bom (0,0). The argument against 

the diagonal 2D spin stripes is that the above charge 
peaks have not been observed, while, on one occasion, 
the charge peaks expected for the ID superstructure were 
observed 

Finally, the third and, perhaps, the most direct exper- 
imental observation suggesting the one-dimensionality of 
the stripe superstructure comes from Ref.™ (also sup- 
ported by Ref.''^), where it was found that partial de- 
twinning of a YBCO sample leads to a very strong asym- 
metry of the (tt, tt) peak splitting. This asymmetry is in 
quantitative agreement with the expectation that each 
of the twin domains in YBCO has only one kind of one- 
dimensional stripe pattern, and, therefore, partial de- 
twinning should lead to a significant redistribution of in- 
tensities between the four peaks. 

The arguments counter-balancing the above experi- 



mental observations can be the following: 

The first observation, although important, is only in- 
direct and thus cannot substitute for direct observations. 

The interpretation of the second observation is based 
on the assumption that the spin structure is weakly mod- 
ulated by function cos[^(a:: -I- y)] cos[^{x — y)]. As a 
result the charge modulation is also assumed to be small, 
and, therefore, the Landau theory is applied. However, 
if the modulation of the the AF structure is strong, and 
the stripes are, indeed, deep and narrow as assumed 
in Section |nl then the four charge peaks proposed as 
an indicator of the 2D nature of stripe pattern would 
only be a part of a more complicated peak structure, 
and not necessarily the most pronounced one. For 
example, for a more realistic charge profile modulated as 
''-':{x + y)] sm^[^{x-y)] 



siu^l^ix- 



y)] 



cos 



the positions of the strongest satellite peaks coincide 
with those expected for the ID picture. 

The effect of the 45 degree rotation of some modulation 
peaks can, nevertheless, be relevant in another context. 
It will be shown in Section IVIII CI that this effect can 
lead to the checkerboard pattern observed by STM in 
the local density of states of Bi2213i2ii2i2a. 

Concerning the third argument against the 2D diago- 
nal stripe picture, the experimental data, as they stand, 
appear quite convincing. Yet, this kind of evidence has 
been limited so far only to YBCO. The significant in- 
ternal anisotropy of the lattice structure in YBCO can, 
in principle, induce an anisotropic INS response even 
of a truly two-dimensional stripe superstructure. Thus 
the final resolution of the dilemma between 2D struc- 
ture with anisotropic properties and a ID structure can- 
not be relied on this experiment alone. [Later note: A 
very recent INS study^ has shown that the splitting of 
the (tt, tt) peak measured on detwinned YBCO crystals is 
anisotropic but clearly two-dimensional.] 



IV. DYNAMIC PROPERTIES OF THE STRIPE 
STRUCTURE 

Now I address the question, to what extent the stripe 
picture described in Section^]can possess dynamic prop- 
erties. 

There are two possibilities for the time dependence of 
that structure. 

The first possibility is that the stripe boundaries of 
AF domains can fluctuate and then drift away. One cir- 
cumstance that favors the fluctuations of stripes with or 
without average drift is that, in the approximation of the 
nearest neighbor exchange and without holes placed in- 
side the stripes, the spins located on the diagonal bound- 
aries of AF domains have zero exchange energy and thus 
are free to flip (see FigQJa)). This would not the be 
case if the AF domain boundaries were oriented along 
the principal lattice directions. 

With holes inside the stripes, it is even easier for the 
stripes to fluctuate locally, but it is more difficult to drift 
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on average. The Coulomb repulsion between different 
stripe elements, the topology of the 2D stripe superstruc- 
ture and pinning on impurities and structural disorder 
should inhibit the average drift of the stripe pattern. 

Thus, in the following, I will assume that, on the 
timescales relevant to the physics of superconductivity, 
the stripe superstructure does not drift, even though in- 
dividual stripe elements may exhibit strong transverse 
fluctuations. 

The second possibility for the time dependence of the 
stripe superstructure is that the AF order parameter of 
a given AF domain can fluctuate and then, perhaps, ex- 
hibit "rotational diffusion". The relative spin orienta- 
tion of neighboring AF domains is flxed not by direct ex- 
change interaction but by the overall energy balance of 
the entire structure, which does not allow holes to leave 
the stripes. Spins of a single AF domain cannot simply 
flip by 180 degrees, because such a flip would dissolve the 
domain boundaries, which, in turn, would contradict to 
the assumption that the high energy balance favors the 
stripe structure. It can, however, happen that the AF or- 
der parameter of a given domain fluctuates slightly away 
from the the anti-alignment orientation with respect to 
the neighboring domains, and then either this fluctuation 
is damped back to the initial orientation, or, on the con- 
trary, the neighboring supercells adjust the orientation 
of their spins, and, in this way, the average orientation 
of the spin order diffuses away for relatively large regions 
of the sample. The range of the spin orientations swept 
by the rotational diffusion will, eventually, depend on the 
relative strength of the local anisotropy interaction with 
respect to the the spin fluctuations. 

As a consequence of the above kind of rotational dif- 
fusion, the elastic response from the modulated AF spin 
structure should be either unobservable or strongly sup- 
pressed. At the same time, I assume that this dif- 
fusion is slow enough and thus can be neglected in a 
theoretical model of superconductivity. (To date, the 
elastic response of presumed stripe superstructure has 
been detected from the superconducting samples of Nd- 
dopeS'^^, Zn-doped^S*^and "pure" LSCO^ and from 
Lai.875Bao.i25-2;Sr2;CuOzf^.) 



V. POTENTIAL PROFILE AND THE TWO 
KINDS OF HOLE STATES 

The potential proflle for the hole excitations in the 
background of 2D stripe superstructure is sketched in 
Fig. 121 This profile consists of two main components: 
the network of narrow potential wells running along the 
stripes and the (repulsive) Coulomb potential created 
around them. 

It is clear from the above picture, that there are two 
kinds of hole states: those localized mainly inside the 
stripes — I will call them b-states — and those local- 
ized mainly in the shallow potential wells inside the AF 
domains — I will call them a-states. 




FIG. 3; (a) Sketch of potential landscape within one stripe 
supercell. (b) Scheme of a stripe supercell. (c) Sketch of 
potential landscape along the thick line drawn across plot 
(a). Energies Ek, Em and Ep correspond to points K, M 
and P marked in figure (b). Solid horizontal lines indicate 
quantum levels inside the respective potential wells. Dashed 
line indicates the position of chemical potential fi. 



It is, particularly easy to discuss the situation, when 
(A) the in-stripe potential wells are deep enough, so 
that, at low temperatures, almost all holes stay inside 
the stripes, and (B) the AF domains are large enough, 
so that the energy levels of both a-states and b-states 
are spaced closely enough, as shown in Fig. |3Ic). In this 
case, the chemical potential fj, of holes is approximately 
equal to the highest energy of occupied b-state at zero 
temperature. This energy should be below Ek (the min- 
imal energy of a-states indicated in Fig. |3Ic)), which is 
a formal restatement of the assumption that the stripes 
are deep. 

Such a situation may well correspond to strongly un- 
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derdoped (but superconducting) cuprates, in which case, 
the difference — fJ- can be identified with the pseu- 
dogap defined as the leading-edge midpoint in angle- 
resolved photoemission spectroscopy (ARPES)^». The 
phenomenology of the pseudogap can then be interpreted 
as follows. At energies between /i and Ek — /i, the den- 
sity of states is low, because only b-states contribute to 
this energy range. Moreover, since b-states are extended 
along the diagonal stripes, the momentum of photoelec- 
trons emitted from b-states should be preferentially ori- 
ented along the lattice diagonals. Above the pseudogap, 
the density of states gradually increases due to the con- 
tribution from a-states. The characteristic energy scale 
of this increase can be estimated as 0.1 — 0.2 eV (a typical 
height of the Coulomb potential barrier Em — Ek shown 
in Fig. He)). 

If the assumptions (A) and (B) hold at higher dop- 
ing concentrations, then one can envision some critical 
concentration (perhaps, not far from the optimal dop- 
ing), at which the pseudogap measured as Ek — fJ. be- 
comes equal to zero. However, the other measure of the 
pseudogap, the gradual increase of the density of states 
above Ek = /i, will not disappear as long as stripes re- 
main stable and, therefore, generate locally inhomoge- 
neous Coulomb potential landscape. 

One should note, however, that the validity of assump- 
tion (B) above is, particularly questionable at the dop- 
ing concentrations, corresponding to the physically in- 
teresting sizes of AF domains of the order I = 4ao-\/2 « 
23 A(for flo = 4 A). If, for an estimate of the level spacing 
of both a- and b-states, one takes the spacing between 
the lowest levels of a free particle having effective mass 
771* — bnie in a box of size I, then one obtains 40 meV. 
(Here 777e is the bare electron mass.) The number 40 meV 
is of the order of the experimental values of the pseudo- 
gap and also notably larger than a typical critical tem- 
perature (« 7 meV). 

If the level spacing is, indeed, as large as estimated 
above, then it is likely, that there are no states (or very 
few due to disorder) in the energy range between the 
highest occupied b-state and the lowest a-state. In such 
a case, the position of the chemical potential within the 
above energy window becomes uncertain. However, the 
pseudogap can still be defined as Eq — where Sa is the 
lowest energy of a-states, and Sb the energy of a b-state 
closest to the chemical potential. 

Formally speaking, the latter definition extends to a 
somewhat counterintuitive situation characterized by in- 
equality Ea < Sb- While, in the underdoped cuprates 
the strong expectation is that > Eb for both large 
and small values of the level spacing, the analysis of ex- 
periments in Section IVIII Dl will suggest that inequality 
Sa < Sb may characterize high- Tc cuprates having doping 
concentrations above Xdc ~ 0.19. The condition Sa < Sb 
contradicts to the assumption (A) made above. However, 
as long as there are only a few a-states with energies 
smaller than ef,, a significant fraction of holes will still 



be localized inside the stripes and generate the Coulomb 
potential required for the validity of the sketch shown in 
Figs.|la,c). 

One should also note that the condition Sa < Sb is 
detrimental to stripe stability, because it implies that 
holes penetrate inside the AF domains. However, the 
emergence or disappearance of stripes is not just the 
subject of one-particle considerations. This process is 
governed by the balance of global energy, which, among 
others, includes the contributions from lattice strain and 
quantum AF fluctuations. It is, therefore, not necessary, 
and, in fact, unlikely, that stripes become unstable pre- 
cisely when Sa = Furthermore, in the present work, 
the condition Sa < £& will only imply one filled a-state per 
AF domain, which the stripe superstructure may, indeed, 
sustain. 

Now, I motivate the form of the model Hamiltonian, 
which will be introduced in the next Section. 

Describing b-states, I assume that they are fermionic 
states carrying charge e and belonging to one stripe el- 
ement. It will only be important for the model that b- 
states can form pairs with total spin 0. Whether or not 
they carry spin 1/2 is not of primary importance. 

An a-state is the state of one hole injected into a fi- 
nite AF domain. It is not important for the model to 
know exactly the orbital and the spin wave functions of 
a-states. It is, however, important to note that (i) the 
spin wave function of an a-state should be fixed by the AF 
background, which is assumed to be static on the time 
scale of interest, and (ii) the AF order parameter has op- 
posite sign for two neighboring AF domains. Therefore, 
if two orbitally equivalent a-states from neighboring AF 
domains form a pair, then the total spin of that pair is 
equal to zero. 

It should be mentioned here, that the analogs of a- 
and b-states have been identified in the numerical study 
of spin polarons in the stripe background'^'^. 

As discussed in Section IIVI the diagonal orientation 
of stripes predisposes them to strong transverse fiuctu- 
ations. These fluctuations should, in turn, strongly in- 
teract with both a- and b-states. The effect of this in- 
teraction is then two-fold. On the one hand, the stripe 
fluctuations couple to holes both elastically and inelas- 
tically, and thus suppress the hole transport across the 
stripe superstructure. On the other hand, they can ef- 
flciently mediate the interaction between different hole 
states. 

The first effect justifies the following "center-of-mass" 
selection rule: The model Hamiltonian can have tran- 
sition elements only between quantum states having the 
same center-of-mass coordinate. 

This selection rule, first of all, eliminates the direct 
hopping terms between a- and b-states belonging to dif- 
ferent AF domains or stripe elements. Since this is quite 
a radical assumption, here I list several additional factors, 
which contribute to the suppression of hopping in addi- 
tion to the stripe fluctuations. These factors are: (i) mis- 
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match of AF backgrounds between two neighboring AF 
domains; (ii) Coulomb potential barriers between neigh- 
boring AF domains and between intersecting stripe ele- 
ments; (iii) disorder in the stripe superstructure (which is 
not present in Fig. |3| but should be present in a real sys- 
tem). For the subsequent treatment, it is not important, 
whether a- and b-states are rigorously localized. All the 
above factors together should only ensure that the hop- 
ping terms are small in comparison with the interaction 
term discussed below. 

In the absence of hopping, every a-state and every b- 
state are to be characterized by "on-site" energies Sa and 
Eb respectively. 

I assume that diagonal interactions involving a- and 
b-states, and also non-diagonal interactions between a- 
states inside the same AF domain and between b-states 
inside the same stripe element can be satisfactorily taken 
into account by the renormalization of energies ea and et- 

Considering the alternatives for non-diagonal interac- 
tion terms, I limit the model choices only to the terms 
of the fourth order with respect to the fermionic creation 
and annihilation operators. Most of these terms fail to 
qualify under the "center-of-mass" selection rule. Among 
the few terms, which qualify, the only one, which will be 
included in the model, corresponds to the transition of 
two holes occupying two a-states on the opposite sides of 
a given stripe element into two b-states inside that stripe 
element and vice versa (see Fig.(0Jb)). One can check 
that the center of mass of two "initial" a-states coincides 
with the center of the stripe element between them, and 
thus coincides with the center of mass of two "final" b- 
states. 

The above transition can be efficiently mediated by the 
fluctuating stripe element itself. The relevant mechanism 
would involve two steps. Step 1: A hole hops between two 
adjacent a- and b-states. Since the two states have differ- 
ent centers of mass, the transition between them should 
be accompanied by a virtual excitation of the transverse 
oscillation mode of the stripe element "housing" the b- 
state. Step 2: The above oscillation mode is absorbed in 
the course of the symmetric transition of a second hole, 
which involves an a-state from the other side of the same 
stripe element and, therefore, restores the center-of-mass 
position. In order to appreciate this mechanism, one can 
look at Fig.|31 and imagine that, in the course of a trans- 
verse fluctuation, one of the two potential wells, which 
represent stripes, shifts towards the center of the figure. 
As a result of this shift, the wave function of an a-state in 
the center of the figure strongly overlaps with the wave 
function of a b-state inside the shifted stripe. Such a 
strong overlap constitutes a precondition for a large value 
of the interaction term selected for the model. 

The above term, which can be schematically described 
as "aa^bb," is sufficient to achieve the primary goal of 
the present work, which is to find at least one plausible 
channel for the superconductivity in the presence of the 
2D arrangement of diagonal stripes. There exist, how- 
ever, a few other fourth order non-diagonal terms having 



form "ab^ab", "aa^aa," or "bb<-^-bb," which qualify 
under the center-of-mass selection rule. At the moment, 
I rank these terms as less important, because either they 
involve states, which are too far separated from each 
other, or they imply no charge fiow between different 
components of the stripe superstructure. Nevertheless, 
the effect of including additional non-diagonal (and also 
diagonal) terms in the Hamiltonian would merit further 
study. 



VI. SUPERCONDUCTIVITY MODEL 
A. Hamiltonian 

In this Section, I introduce a model which is minimally 
sufficient to capture all the elements of the qualitative 
description given in Section Ivl 

The model is limited to the two-dimensional superlat- 
tice of a- and b-states, which is shown in Fig.^ This su- 
per lattice is divided into diamond-shaped supercells, like 
the one shown in Fig. Ol^b) . The supercells are labelled 
by single index i (or j). Each of them is characterized by 
the AF index rji (defined in Section ^] but now with sub- 
script i). I will use the terms "even supercell" and "odd 
supercell" to refer to the supercells having rji ~ 1 and 
rji — —1 respectively. The total (macroscopic) number of 
supercells in the system is denoted by variable N. I will 
further assume that the system has rectangular form of 
dimensions and Ly along the x- and the y-directions 
respectively. 

In the present model, there exists only one a-state in- 
side each supercell described by hole annihilation and 
creation operators Oj and af respectively. Each a-state 
is characterized by the on-site energy Sa- In order to 
implement the observation made in the previous Section 
that orbitally equivalent a-states from the neighboring 
supercells have opposite spins, the spins of a-states al- 
ternate together with the AF index 77^, i.e. the spin of 
ith a-state is equal to ^rji. Operators do not need an 
additional spin index, because, in this model, the spin of 
an a-state is not an independent quantum number but 
instead fixed by the lattice index i. 

The next assumption is that, inside each stripe element 
separating the ith and the jth neighboring supercells, 
there exist only two b-states both characterized by the 
same "on-site" energy £{, but having different spins (+1/2 
or -1/2). The hole annihilation and creation operators for 
these states are and bfj ^ respectively. Here index a 
represents the spin of a b-state and assumes values -I- and 
— . (The assumption of a well defined spin is made just for 
the concreteness of the model. In reality, it can, indeed, 
be spin but also any other kind of quantum number.) 

In total, this model contains one a-state and four b- 
states per one supercell (two b-states per stripe element 
times two stripe elements per supercell). 

The on-site energies £a (the same for all a-states) and 
Eb (the same for all b-states) are both measured from 
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(a) 




(b) 



+ ^ 



a: 



b 



where g is the interaction constant. Here and below the 
notation in the sum subscript imphes that the sum 
over index j extends only over the nearest neighbors of 
the ith supercell. Expression = 1 in the sum super- 
script means that the summation over index i includes 
only even supercells (marked by pluses in Fig^fb)). 

Hamiltonian 10 can be described as an exotic 2-band 
model, where the non-interacting states are localized but 
still have the same on-site energy, and the interaction 
includes only inter-band coupling. In such a model, the 
variational SC ground state exists independently of the 
sign of the coupling constant g. In the calculations, I 
will, therefore, assume that g > 0. 

Two b-states per stripe element and one a-state per 
AF domain, represent the minimal configuration required 
for implementing the interaction term in the Hamilto- 
nian Q. At the same time, it should be noted here 
that the ratio "one hole per 2ao\/2" extracted in Sec- 
tion ^ from the scaling of INS data also corresponds to 
2 holes per stripe element in the most interesting case of 
I fa 4ao-\/2. Furthermore, if the level spacing inside AF 
domains is, indeed, as large as estimated in Section Ivl 
then keeping only one a-state per AF domain also con- 
stitutes a meaningful approximation for the description 
of low-energy properties of the system. 



B. Classification of the model regimes 



+ 



A 



The present work is mostly limited to the analytical 
results describing the following regimes: 

Case lA: Sa > 0, Sb = 0; 



FIG. 4: (color online) (a) Model quantum states: one a-state 
with energy Sa inside every AF domain, and two degenerate 
states with energy et inside every stripe element. This picture 
is to be compared with Fig.^^c). (b) Two-dimensional scheme 
of a-states and b-states. Each circle represents the center 
of an a-state, while each ellipsoid extended along the stripe 
boundaries represents the location of two b-states. Signs 
or " inside the circles indicate the sign of AF index rji. 
Also shown: two operators of a-states and two operators of 
b-states in their respective spatial domains. The transition 
between these two pairs of a-states and b-states represents a 
typical term in the interaction part of Hamiltonian Q . 

the chemical potential. More detailed assumptions about 
their values will be made later. 
Finally, the model Hamiltonian is: 

+ 9^2 i^t+^t-'^^a, + h. c), (7) 



Case IB: Sa <0,Sb~ 0; 
Case HA: = 0, ef, > 0; 
Case HE: Sa = 0, < 0; 
Critical case: Sa ^ £b = 0. 

The diagrams of energy levels representing Cases lA, 
IB, HA and IIB are sketched in Fig. |S1 The reasons for 
distinguishing those special regimes from the general case 
{Sa 7^ 0, £{, 7^ 0) are the following: (i) As shown later in 
subsection lVI Dl for the fixed difference Eq — e^, the situ- 
ations, when cither Sa or Sb coincides with the chemical 
potential (i.e. equal to zero) correspond to sharp minima 
in the ground state energy, (ii) The condition Sa — or 
Eb = leads to a significant simplification of the model 
calculations. 

Since most observables characterizing Cases lA and IB 
are identical, these two Cases will be referred to as Case I, 
whenever the difference between lA and IB is not impor- 
tant. Similarly, "Case H" will refer simultaneously to 
Cases HA and IIB. 
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IIB 




FIG. 5: Sketches representing Cases lA, IB, II A and IIB dis- 
cussed in the text. 



C. Mean- field solution 

In this Subsection I shall proceed with finding the vari- 
ational ground state using the method of Bogoliubov 
transformation. The variational procedure will consist 
of (i) making Bogoliubov transformation for b-states in 
real space; (ii) truncating the Hamiltonian Q by leaving 
only the diagonal terms with respect to the new Bogoli- 
ubov quasiparticles and then averaging those terms; (iii) 
making the Fourier transform of a-states; (iv) introduc- 
ing the Bogoliubov transformation of a-states in fc-space; 
and, finally, (v) minimizing the energy with respect to 
the both transformations. Although straightforward, the 
above procedure is somewhat involved. For this reason, 
in Appendix^ I also present an approximate version of 
the same mean-field solution. The content of Appendix IXI 
also reveals a number of interesting facts about the ro- 
bustness of the full solution. 

The first step of the variational scheme consists of the 
following Bogoliubov transformation: 



= sSy_ - we'^^'^B. 



(8) 
(9) 



where i?^,^ are the annihilation operators of the Bogoli- 
ubov quasiparticles, s and w are the (real) transformation 
coefficients satisfying the normalization constraint 



1, 



(10) 



and ifij are the transformation phases chosen to be the 
same for all translationally equivalent stripe elements. 

There exist four translationally non-equivalent types of 
stripe elements. Each type corresponds to one of the four 
possible orientations of vector rj—ri, where and rj are 
the positions of the centers of two neighboring supercells. 
One of these two supercells is always even, while the other 
one is always odd. I will use the convention assigning 
to an even supercell, and r, to an odd one. The four 



possible realizations of vector Vj — are: 

R. = ^(-i.i); 

R. = -ij,l.-l). 

Phases ipij can now be presented as 
tpij = ip{rj - Ti). 



(11) 
(12) 
(13) 
(14) 

(15) 



They can have, at most, four different values ipi — (/^(Ri), 
If 2 = <i2(R2), — 'Pi'^s) and = </j(R4). I label 
these four phases by index a and refer to them using the 
notation 



fa 



(^(Ra). 



(16) 



The physical explanation, why the four phases (fa should 
be tracked in the variational solution, is given in the end 
of Appendix ^ 

Substituting transformation JS] EJl into the Hamilto- 
nian Q, and then averaging the result with respect to 
B-operators, I obtain 



Ha = 'ieb N [s^HB + w^(l - ub)] +£aY^ a+a. 



-\-gsw[l — 271b) ^ 



')i=i 



h. c. 



where 



Ub 



exp(^) +1' 



,(17) 



(18) 



Here, ns and Eb are, respectively, the occupation num- 
ber and the energy of a B-quasiparticles; and T is the 
temperature measured in energy units. Index j(i, Rq.) in 
Ea. (|17|) labels the nearest neighbor of ith supercell such 
that Yj — Yi = Rq . 

At this point, it is convenient to replace supercell in- 
dices i and j in Ea. (|17l) by the set of the radius- vectors 
and also to separate explicitly the summations over even 
and odd supercells. This gives 



TLa = 4e&7V [s^uB + w^(l - ub)] 
+ea X! a+(re)a(re) + a+(ro)a(ro) 



(19) 



re 



-gsw{\ — 2nB) 



"■o 



J2 a{Ye)a{Ye + R„) + h.c. 
re 



where even supercells are characterized by the discrete 
set of radius- vectors {re} and odd supercells by the com- 
plementary discrete set {ro}. Note: any vector of the 
form re + Rq belongs to the "odd" subset. 
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Now I introduce the Fourier transform separately for 
even and odd supercells: 



re 

ro 



-ik: 



ro 



(20) 
(21) 



The two Fourier transforms, although involve different 
parts of real space, still performed with the same set of 
fc- vectors, because the even and the odd subsets have the 
same periodicity. The projections kx and ky of the k- 
vectors, change in discrete steps ^ and respectively. 
They fall in the interval — § < fc^ , Aij^ < , where d is the 
period of the sublattice of even (or odd) supercells equal 
to l\/2. The total number of A:- vectors is 



Nk = LxLy/{d?) = N/2. 



(22) 



After transformation H20I21|I . the Hamiltonian H19|) can 
be written as 



^ebN[s^nB + w'^{l-nB)\ (23) 
+ea ^ at (k)ae (k) + Sa'^a'^ (k)ao (k) 

k k 

+ gsw{l — 2nB) 

x^[ae(k)ao(-k)F(k)+h.c.], 



where 



-ikRa 



nk) ^Ea^-'^' 

= 2 exp [— 
+2 exp [-i ^22±^ 



i cos 



[kRi 
cos [kR2 



V2—Vi 1 
2 J ■ 



(24) 



Given the form of the interaction term in the Hamiltonian 
(|23|l . the choice of canonical transformation for a-states 
is clear: 

ae(k) = u{k)Ae{k)+v{k)e'^-^''^A+{-k), (25) 
ao(-k) = u{k)Ao{-k)-v{k)e'*-^'''>A+{^), (26) 

where Ae(k) and Ao{k) are annihilation operators of the 
new Bogoliubov quasiparticles; 0a (k) is the phase of this 
transformation; and u{k) and v{k) are the real numbers 
obeying the following normalization condition: 



u^{k)+v^{k) = 1. 



(27) 



An important conceptual detail to be noted here is 
that transformation (|25l I26|l will eventually lead a coher- 
ent one-particle dispersion of A-quasiparticles in fc-space. 
This kind of fc-space coherence emerges only in the SC 
phase and appears to be "protected" by the Fermi statis- 
tics (see the end of Appendix 



Substitution of transformation H25II26|I into the Hamil- 
tonian H19|) results in the following expression for the en- 
ergy of the system: 



E = 4:etN[s'^nB + w'^{l~nB)] 



(28) 



+2ea 2^ {u\k)nAik) + v\k)[l - ^^(k)]} 

k 



+ 2gsw{l - 2nB) 
u{k)v{k){2nA{k) - 



l)|V^(k)|cOs[0y(k)+0a(k)], 



where |V^(k)| and (/)v(k) are the absolute value and the 
phase of the complex- valued function l|24|l . and 



nA{k) 



exp 



(^) 



1 



(29) 



Here ^^^(k) and £A(k) are, respectively, the occupation 
number and the energy of a Bogoliubov quasiparticle cre- 
ated by operator (k) or Aq (k) . 

It is immediately obvious that the minimiza- 
tion of the above expression requires the last 
term to have maximally negative value. All the 
sign conventions used below will be such that 
gsw{l - 2nB)u{k)v{k){2nA{k) - 1) < 0. Therefore, the 
maximally negative value will be reached, when 



COs[0y(k) -f </.a(k)] - 1, 



(30) 



which means that (up to an integer number of 27r's) 

</)a(k) = -0y(k). (31) 



For the reasons discussed in subsection IVI Bl the rest 
of the calculations will be mostly limited to Cases lA, 
IB, HA and IIB. Formulas for the critical case will be 
summarized in subsection IVI Gl For the general case, 
only the condensation energy will be obtained (in sub- 



Case lA: Eb = 0,ea > 0. 

A natural sign convention in this case is: eyi(k) > 
and es < 0, i.e., at T = 0, n^(k) = and hb — I- 

Given the above convention, the minimization of en- 
ergy gives 



u(k) 
v{k) 



\ 



1 1 




\ 




where 



Q(k) 

r(k) 



1 

' 71' 
1 

2eai2nA(k) - 1); 
g(l-2nB)(2n^(k)-l)|y(k)|; 



(32) 

(33) 

(34) 
(35) 



(36) 
(37) 
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By varying E with respect to uaO^) and ns, one can 
now obtain the quasiparticle energies: 



1 



-gH2nB~ir\VikW 



(38) 



.9^(2nB-l) ^ (l-2ru(k)|nk)P 
= ^7(k^ ' ^''^ 



8A^ 



and then express the total energy of the system as 

E = ^Y.^{l-2nA{k))eA{k)-ea]. (40) 

k 

Here and everywhere, the summation over k can be re- 
placed by integration according to the following rule 



1 — 

iV ^ ^ 87r2 



(41) 



Energy 140(1 is the function of |l^(k)| (via eyi(k) and 
nA(k) = n^[e^(k)] ). In turn, |V^(k)| is the function 
of four phases ((16|l (via Ea. H24|l 'l. Therefore, energy 140|l 
should further be minimized with respect to the values 
of those phases. 

As shown in Appendix^ such a minimization imposes 
only one constraint: 



= — h 7rn, 

2 2 ' 



(42) 



where n is an integer number. Among three other inde- 
pendent combinations of phases (pi, Lp2, and Lp4, one 
should remain free as a consequence of the global gauge 
invariance, while two remaining combinations should, 
in principle, be fixed, but not in the framework of the 
present model. 

The main thermodynamic and transport properties of 
the model are independent of the choice of phases ipi, 
^2, f3 and ifi as long as this choice is consistent with 
Ea. (|42ll . This situation is somewhat similar to that of 
superfluid "^He, where the interaction causing the super- 
fluid transition does not fix the values of all variables 
characterizing the order parameter—. In '^He, the re- 
maining freedom is eliminated by magnetic dipolar inter- 
action between nuclei, and by other small interactions. 
In the present case, the same role can be played, e.g., by 
pair hopping between a-states or b-states belonging to 
different supercells. In this work, the issue of the "phase 
freedom" is not resolved. It is, however, possible to spec- 
ulate that the additional terms will lead to sufficiently 
symmetric selection of phases, such that 

V{k) = 2 {cos [kRi] - i cos [kRz]} , (43) 



or 



l/(k) = 2{sin[kRi] -i sin[kR2]}, (44) 



V{k) = 2 {cos [kRi] - i sin [kRs]} 



(45) 



The first choice corresponds to ipi = tpz = 0, = y^i ~ 
7t/2; the second one to (pi — —tt/2, (p2 — 0, (p3 — tt/2, 
and (p4 = tt; and the third one to ipi — 0, ip2 — 0, ips — 0, 
and (/34 = tt; The resulting patterns of phases are shown 
in Fig. ®. 

Having specified the parameters of transformations 
El Uni [2n|l ; one can calculate the temperature dependence 
of various thermodynamic quantities. This requires the 
numerical solution of the system of equations H181 1291 1381 
I39|l . which is not done in the present work. Without the 
full numerical solution only the zero-temperature char- 
acteristics and the SC transition temperature can be 
evaluated. The evaluation of Tc is based on a manipula- 
tion described in Appendix |21 which gives the following 
simple equation: 



9' 


exp (f^^ 


- 1 




exp (f^^ 





(46) 



When g < Ea, the approximate solution of Ea. H46|l is 
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Sea' 



(47) 



Another simple limit is Sa = 0, in which case, Eq. I|46() 
yields 



Tr=^ 



(48) 



In general, however, Eg. lfl^ has to be solved numerically. 

Since the operators of both real holes and real elec- 
trons do not commute with A- and B-operators defined 
by Eas. l|AllA3l) . the tunneling studies of both A- and B- 
quasiparticles (via a contact with normal metal) should 
show the density of states on the both sides of the chem- 
ical potential, i.e. at e = ±e^ and e = zLeb- Moreover, 
as long as et = 0, tunneling into B-states should result 
in the density of states symmetric with respect to the 
chemical potential. As far as A-states are concerned, 
then tunneling into them should show asymmetric den- 
sity of states. This asymmetry is characterized by the 
ratio 



D[eA{k)] _ u^k) _ eA{k)+ea 



D[-eAik)] v^{k) eA{k)~ea 



(49) 



The zero-temperature tunneling spectra of A- and B- 
quasiparticles are shown in Fig.jTfa). The spectrum of B- 
quasiparticles consists of two delta peaks located at ±63 ■ 
A-quasiparticles have a continuous spectrum, which is 
fully gapped with minimal energy Sa- It was obtained 
by first calculating the density of states following from 
Ea. (|38|) as a function of positive energies ea and then 
dividing the weight between positive and negative tun- 
neling energies according to formula H49|) . 
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FIG. 6: (color online) Three examples of particularly symmet- 
ric patterns of phases ipij consistent with the variational SC 
solution. The values of the phases are indicated on the top of 
the corresponding stripe elements. Each pattern is obtained 
by the periodic translation of four circled phases denoted in 
the text as (pa- These phases are constrained by Ea. (l42ll . The 
expressions for V^(k) corresponding to the phase patterns (a), 
(b) and (c) are given, in respective order, by Eas. l43L I44II 
and lH^ . 



-e/g 



FIG. 7: (color online) Examples of the tunneling density of 
states at T = 0: (a) Case lA {sa = 0.5g, et = 0); (b) critical 
case {sa — 0, Eb — 0); (c) Case IIA [sa — 0, et = 0.5g). In 
each case, the spectra contain two vertical delta-peaks rep- 
resenting B-states and located at ±es given by (a) Eg. 1381 : 
(b,c) Eg. 15911 The continuous part in each spectrum represents 
A-states. It is calculated from: (a) Eg. 1381 : (b,c) Eg. 1591 . In 
all three cases, the spectra of A-states have Van-Hove singu- 
larities located at ±eAo and the termination points located 
at ztsAi- The spectrum of A-states in figure (a) also has 
a gap between Sa and —Sa- The asymmetry of the spectra 
is obtained from: (a) Eg. 1491 : (c) Eg. 1641 . Note: the posi- 
tive direction of the horizontal axis corresponds to negative 
hole energies. (This reflects a convention of tunneling spec- 
troscopy.) 



In addition to the gap and the asymmetry, two other 
important features of the spectrum of A-quasiparticles 
are: the Van-Hove singularity and the sharp termination 
point at a higher energy. These two features correspond, 
respectively, to the saddle points and to the maxima of 
|y(k)|. Function |l^(k)| obtained from Ea. (|43|l is shown 
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(a) 



V(k)| 




Altman and Auerbach from the plaquette boson-fermion 
model^^. 

Case IB: et = 0, < 0. 

In this case, if one keeps the same sign convention as in 
Case lA, then u(k) and v(k) given by Eqs.|j221IS3) should 
switch values. All formulas for the quasiparticle ener- 
gies and the tunneling spectrum asymmetry obtained for 
Case lA apply without modification to the present case. 
Equation H4()|) for the critical temperature also applies 
but with the trivial substitution of |£a| instead of Eq. 

The only observable difference between Cases lA and 
IB is the opposite asymmetry with respect to the chem- 
ical potential: in Case lA, the density of A-states is 
greater on the hole side, while, in Case IB, on the electron 
side. 



(b) 



7T 




FIG. 8: (a) Three-dimensional plot and (b) contour plot of 
|l^(k)| corresponding to Eg. 14311 . Each plot shows the first 
Brillouin zone of the stripe superstructure. In the contour 
plot, the saddle points are indicated with "5"', and zeros with 
"O" . 



in Fig.|Sl It has four saddle points at ks — ^(±1,±1). 
With another choice of phases consistent with Ea. (|42(l . 
the fc-space position of the saddle points may change but 
not the value of |y(ks)| = 2. Therefore, according to 
Ea. (|38|l . the density of A-states exhibits a Van Hove sin- 
gularity at 



(50) 



The maximum energy of A-states, which corresponds to 
the tunneling spectrum termination point, can be found 
by substituting the maximum value of |I^(k)| equal to 
2\/2 into Eq.(|2Hl, which gives 



SAi = ^Jel + 2g\2nB-V? 



(51) 



It should be noted here, that the appearance of the 
tunneling spectrum of A-states shown in Fig. Eta) is 
quite similar to that of the fermion spectrum obtained by 



Case IIA: Sa = 0, £& > 0. 

The sign convention in this case is chosen to be oppo- 
site to that of Case lA, namely: eA(k) < and > 0, 
i.e., at T = 0, n.^(k) = 1 and tlb = 0. 



In this case, 



1 

71' 
1 

71' 



1 1 



T2 



W 




where 



Q = 47V£fc(2nB - 1); 
r = .g(l - 2nB)CaN- 

Ca = l5](2n^(k)-l)|y(k)| 



N 

k 

One can then obtain the quasiparticle energies 
<?2 |y(k)| Ca (l-2nB) 



8£i 



4 + Yq9' 



and the total energy 

E = -2N[{l-2nB)eB-eb]. 



(52) 
(53) 

(54) 
(55) 



(56) 
(57) 

(58) 



(59) 
(60) 

(61) 



The choice of phases (jl6|l is still constrained by con- 
dition H42|l (see Appendix 151). Given this constraint, the 
zero temperature value of Ca (obtained numerically) is 



Caa 



l^|I/(k)| = 0.958., 



(62) 
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A manipulation analogous to the one described in Ap- 
pendix[n]gives the following equation for the critical tem- 
perature: 



9' 


exp (f^) 


- 1 


8eb 


exp 





(63) 



The tunneling density of states corresponding to 
Eas. H59l60|) is shown in Fig. (TJc). Contrary to the result 
for Case I, the tunneling density of A-states in Case II is 
symmetric, while the density of B-states is asymmetric. 
This asymmetry is characterized by the ratio 



DjeB) 



(64) 



An important feature of Case II, which is absent in 
Case I, is that the energy spectrum of A-quasiparticles is 
gapless with the linear density of states around the chem- 
ical potential. Indeed, saO<-) given by Ea. lfS^ touches 
zero in an isolated set of non-analytic points correspond- 
ing to |V^(k)| = 0. For the specific choice of V^(k) given 
by Ea. H43(l . the zeros of |l^(k)| are shown in Fig.|Hl They 
are located at ko = f (±1,0) and kg = f (0,±1).) This 
feature is a direct consequence of the phase relation (021 • 
A deviation from that relation would produce a line of ze- 
ros, which implies a non-zero density of states at sa = 0. 

The density of A-states in Case IIA has Van Hove sin- 
gularity and the termination point located, respectively, 
at 



Sao 



4es 



and 



£A1 



V2 



£A0- 
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Case IIB: Sa = 0,et< 0. 

All formulas for the quasiparticle energies and the tun- 
neling spectrum asymmetry obtained for Case IIA apply 
without modification to Case IIB. Equation H63|) for the 
critical temperature only requires the substitution of |et,| 
instead of et- The only difference between Cases IIA and 
IIB is the opposite asymmetries in the tunneling spec- 
tra of B-quasiparticles. In Case IIA, the hole side of the 
B-quasiparticle spectrum has more weight, while in Case 
IIB, the larger weight is on the electron side. 



D. Chemical potential as a variational parameter 

In this subsection, I argue that the situations corre- 
sponding to fi — Eb (Case I) or /X = Sa (Case II) should 
be considered as likely scenarios describing realistic stripe 
systems. 

The constraint on the total number of particles, which 
is usually used to fix fj,, cannot be straightforwardly ap- 
plied to the present model for the following reasons: (i) 



The model quantum states, form a subset of all quantum 
states of a real "striped" system, and, therefore, the ac- 
tual total number of particles cannot be reliably counted, 
(ii) The system can always readjust the periodicity of the 
stripe superstructure, which would change the ratio be- 
tween the number of model states and the number of 
holes doped into Cu02 planes, (iii) The chemical po- 
tential can change within the model pseudogap without 
affecting the total number of particles occupying model 
states (at T = 0). 

It is, therefore, reasonable to treat the chemical po- 
tential as a variational parameter, which is fixed by the 
minimization of the total energy of the real system con- 
sidered as the sum of the contribution from the model 
states and the contribution from environment (unspeci- 
fied here). 

The model contribution to the total energy as a func- 
tion of the chemical potential can be obtained by solving 
the model in the general case: Sa 7^ 0, £{, 7^ 0. The de- 
scription of general case is as straightforward as that of 
Cases I and II. However, the minimization routine pro- 
duces an integral equation, which couples the values of 
M(k), w(k), s and w, and which has to be solved numeri- 
cally. 

Figure IHl shows three representative curves for the evo- 
lution of the SC ground state energy as a function of the 
chemical potential. Each curve was obtained numerically 
for the fixed values of Sa, £b and g indicated in the cap- 
tion. In order to allow for the variation of /i, the reference 
point for one-particle energies was shifted (in this part 
only) from fi — to some arbitrary value £f, = 0. The 
SC ground state energy was measured from the energy 
of the normal state. The absolute value of thus defined 
quantity is conventionally called condensation energy. 

In a generic situation (two lower curves in Fig. |^ , the 
ground state energy has two minima: the deeper one — 
corresponding to Case I, and the shallower one — corre- 
sponding to Case II. For a fixed value of 5, the SC transi- 
tion is possible only, when the value of ^ lies close enough 
to either Sa or Sb- Dependently on the ratio \ea — £b\/9, 
the values of /i compatible with superconductivity fill ei- 
ther one finite interval including both Sa and Sb or two 
disconnected finite intervals around Sa and £{,. The top 
curve in Fig. |3 illustrates a non-generic situation, when 
the two minima coincide. The resulting single minimum 
then corresponds to the critical case. 

The main purpose of Fig. I^is to illustrate, that, not 
only Cases I and II correspond to local minima in the 
/i-dependence of the ground state energy, but also that 
these minima have the form of cusps. This implies that, 
if the environment contribution to the total energy of the 
model-plus-cnvironment system is a smooth function of 
fi, then the total energy should have the cusp minima 
due to the model contribution at exactly the same val- 
ues of fi as those obtained without environment. If only 
the model contribution were taken into account, then it 
would follow from Fig. |51 that the minimum correspond- 
ing to Case I is the global one. However, the environ- 
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FIG. 9: Condensation energy (with the negative sign) as a 
function of the chemical potential fi. Each of the three curves 
was calculated on the basis of Hamiltonian (|7J by fixing Sa, 
Eb and g and then varying /x. In the calculations, g and et 
were equal, respectively, to 1 and in all three cases, while Sa 
admitted the following values (top to bottom): 0, 0.6 and 1.1. 
The solid/dashed lines indicate the SC/non-SC ground state. 
The deeper minima of the lower two curves correspond to Case 
I, and the shallower ones to Case II. The single minimum of 
the top curve corresponds to the critical case. The vertical 
coordinates of the middle and the bottom curves are shifted 
by -0.6 and -1 respectively. 



ment contribution can change the relative values of ener- 
gies, corresponding to the two minima and thus shift the 
global minimum to that of Case II. 

In principle, it is also possible, that environment intro- 
duces an extra energy minimum in addition to the two 
cusps described above, and, moreover, that additional 
minimum is the global one. 

The resulting possibilities can be summarized as fol- 
lows: Cases I and II of the present model describe two 
different SC states corresponding to two different posi- 
tions of the chemical potential. In the "best case" sce- 
nario (which I also consider more likely) , one of these two 
states corresponds to the global energy minimum of the 
real system, while the other one represents a metastable 



state. In the "worst case" scenario, both SC states are 
metastable, but the true SC state may still be describable 
by the present model with Sa ^ and Sb ^ 0. 

It is further possible that the position of the chemical 
potential near the sample surface is different from that 
in the bulk. Therefore, the SC state, which is metastable 
in the bulk may become stable near the surface and vice 
versa. In principle, it is also possible that the system 
phase separates and forms domains describable either by 
Case I or by Case II. 

In Section IVIII El I will try to discriminate between 
the SC states corresponding to Cases I and II by mak- 
ing comparison between the model predictions and the 
experiments. 

E. Anomalous correlation functions 

Bogoliubov transformations (|8I9I25I26|I imply that, be- 
low Tc, the following anomalous correlation functions 
have non-zero values: 

*a(k) = (ae(k)ao(-k)) = «(k)«(k)e*»(''' [2nA(k) - 1], 

(67) 

^'6(ry,rp„) = {bij-bpn.+) = swe"^'^" (l-2nB)(5(rij-rp„), 

(68) 

where represents the positions of the centers of stripe 
elements, and 6{. . .) is defined as Kronneker delta on the 
discrete superlattice, i.e. it is equal to 1, when its ar- 
gument is zero, and otherwise. The two correlation 
functions Ht)7|) and Ht)8|) are the two components of the 
SC order parameter corresponding to a- and b-states re- 
spectively. 

In the real space, the first component can be defined 

as 

^a{r^,rJ) = {a{r,)a{rj)). (69) 

It has non-zero values only when its two arguments cor- 
respond to the supercells of different kind (i.e. even and 
odd). The formal structure of \l/a(ri, r^) can be expressed 
as follows: 

^a{re,ro) = ^ ^ *,(k)e*('-e-"-o), (70) 

k 

^'a(ro,re) - -*a(re,ro), (71) 



vl/,(re,r^) = *^(ro,r^) = 0. (72) 

where *a(k) is given by Note: Eq.(|7Tl) is the 

consequence of the fermionic anticommutation rule. 

The coherence length of the order parameter ^'^(ri, r^) 
should be inversely proportional to the one characterizing 
y(k) in fc-space. The examination of Ea. H24l) reveals that 
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the characteristic scale of V^(k) is 7r/|Ri| = tt/I. There- 
fore, the coherence length associated with '^ai^i^'^j) can 
be estimated as the supercell size I. It is likely, that on a 
longer scale ^a(ri,rj) exhibits an oscillatory power law 
decay with the period of oscillations being of the order 
oil. 

The coherence length associated with ^j, is equal to 
zero, which means that only b-states belonging to the 
same stripe element form coherent pairs. 

Two useful quantities, which will later be required in 
the calculation of supercurrent are: 



it can be described as having orientational s-wave sym- 
metry with additional strong phase dependence on the 
center-of-mass position of the paired holes. 

The symmetry of ^'a(i"i, is more different from con- 
ventional analogs. It includes the strong dependence on 
(fij, which traces the phase dependence of 'I'f,(ry , rp„), 
but, in addition, \E'a(ri,rj) also exhibits a sign change 
under the translation by one period of the stripe super- 
structure. This sign change reflects the switching be- 
tween even/odd and odd/even order in Ea. (|71|l . 



(73) 



F. Supercurrent and the penetration depth 



and 



(74) 



where the superscript "n.n." indicates that indices i and 
j represent the nearest neighbors. 

In Case I, the explicit expression for can be ob- 

tained by substituting the values of s and w given by 
Eas. H34l35|) into Ea. (|^ for = rp„, which gives 



(75) 



One can then obtain aUj) by making use of the fact that 



int 



AgN 



(76) 



where E\^^ is the interaction part of the energy (i.e. 
^'int = 2gsw{l — 2rt_B) ...). After i^jnt is evaluated 
with the help of Eas. H32l35|l . one can use Eqs.ltTSlEHl) to 

obtain: 



g{l - 2nB)me'^^^ V- [1 - 2riA(k)]|y(k)|2 



8iV 



k 



£A(k) 



(77) 

The role of index rji in Ea. (|77|l is to supply factor 1 or 
— 1 dependently on whether the first index of a{ij) cor- 
responds to an even or an odd supercell. 

In Case II, the expressions analogous to H77I75|I are 



a-iij) 



gCge'"^'^ {1 - 2nB) 
Sen 



(78) 



(79) 



The essential elements of the symmetry of the SC order 
parameter (|75I77|I or, alternatively, I|78I79|I are captured 
in Fig.|H| This symmetry is characterized by the pattern 
of phases ipij indicated on the top of each stripe element, 
and, in addition, by the pattern of index rji. 

The order parameter 4'5(rjj,rp„) has no dependence 
on the relative orientation of r,, and r„„. Therefore, 



The superconducting properties of the present model 
are unusual, because the superconducting phase stiffness 
comes solely from the interaction term of the Hamil- 
tonian Q. This term induces the fundamental "inter- 
nal supercurrent" associated with the transfer of particle 
density from a-states to b-states and vice versa. The 
translational supercurrent then appears as a gradient of 
the internal one. 

The operator expression for the internal current from 
ith a-state to the surrounding b-states can be obtained 
by considering the time derivative of the particle density 
operator: 



J, 



ab{i) 



]9_ 

' h 



h.c. 



(80) 



Here and everywhere in this subsection, index i corre- 
sponds to an even supercell, and index j to an odd one. 

The sum in Eq. H80|l has four terms — each correspond- 
ing to the transfer of the particle density from the «th AF 
domain into a nearby stripe element labelled by the pair 
of indices ij. Therefore, the operator of translational 
current through the iih supercell (to be denoted as J*) 
can be obtained by assigning the direction to the flow of 
particle density associated with each of the above four 
terms, i.e. 



3{i) 



(81) 



where ny is the unit vector in the direction from the ith 
to the jth supercell. 

The translational current is created, when the proba- 
bility of an a-particle to hop into one of the surrounding 
stripe elements is greater in one direction than in the op- 
posite one. For this reason, the translational current can 
only be carried by a-states. The number of b-particles 
hopping on the both sides of a given stripe element is 
the same for each quantum transition generated by the 
Hamiltonian j?)). 

Now, I show that the phase, which drives the internal 
supercurrent, is 



-<^,(k)-0y(k). 



(82) 
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In equilibrium, according to Ea. H31|l . 4>ab = 0. The 
state with non-zero (j^at can be obtained by modifying 
the Bogohubov transformation for b-states (Eqs. lIHlEl)) 
as foUows: 

= sB,,+ + w;e'('^-+^-)B+_; (83) 
hj- = sB,,- - u;e*('^-+^-)B++, (84) 

while keeping the transformation H25I I26|l for a-states un- 
changed. In this case, the anomalous correlation func- 
tions are, for a-states belonging to the neighboring su- 
percells. 



a{ij)] 



(85) 



and, for b-states belonging to the same stripe element. 



(86) 



where 5'a(ij) and are the equilibrium values of the 

two SC components given, in Case I, by Eas. H75l77(l . and, 
in Case II, by Eas- fTSIT^ . 

The averaging of operator (|8()|l with the subsequent 
substitution of Eas. (|85l8t)|l gives the internal supercur- 
rent: 



{Jab(2)) 



2 g sin0nb 



Z.*bfe)*afe). (87) 



Each term in the above sum is a real number given by 
Ea. (|76|l . (Spatially homogeneous internal supercurrent 
(|87|l can exist, when the particle density oscillates be- 
tween a- and b-states.) 

The translational supercurrent emerges, when (pab be- 
comes position-dependent [to be denoted as (jiabi'^ij)]- In 
this case, the relevant Bogoliubov transformation for b- 
states is 



= sB,. - w e*['^.. (r^)] R+ 



(89) 



The transformation for a-states is still given by Eas. pSI 
I26|l . The averaging of Ea. l81|l under the assumption that 
phases (t>ab{'*^ij) are small and have weak positional depen- 
dence, gives the following expression for the translational 
supercurrent: 



'^*6(y)*a(«i)V(/)a6 



(90) 



Note: according to Ea. (|76|l . the value of the product 
'^b(ij)'^a{ij) is independent of the orientation of the stripe 
element labelled by indices ij. 

From Ea. (|90|l . the supercurrent density can be ob- 
tained as 



Izo 



where 



(91) 



(92) 



zq is the transverse distance per one SC plane, and e the 
charge of electron. As follows from Ea. (|91|l . the uncon- 
ventional feature of the present model is that the super- 
current is induced not by the phase gradients of '^a and 
^'fc separately, but by the gradient of the phase difference 
between and 5*6. Keeping up with convention in the 
literature, I will refer to the phase stiffness as "super- 
fluid density," but the intuitive associations with some 
kind of real density would be misleading in this case. 

In Case I, the substitution of Eqs.(|7S||73 into Ea. l|^ 
gives 

eff^(l-2n^)^ ^ (l-2n^(k))|y(k)p 

- i6A^;.zo ^ T^kj ' 

and, in Case II, the analogous equations 1)781 179|l lead to 
65^(1 - 2nB)Cl 



(94) 



Now I calculate the in-plane penetration depth A of 
magnetic field directed perpendicularly to the SC planes. 

The natural expectation is that the gauge-invariant 
generalization of Ea. (|91|l accounting for the presence of 
the vector potential of electromagnetic field A has form 



[ V(t)ab - ^A 



(95) 



where c is the speed of light. In the present work, I do 
not derive Eq. 1)95(1 but take it as an additional postulate. 

The equivalent of the London limit in the present 
model is A ^ L In this limit, the standard resulli^ for 
the penetration depth, which follows from Eq.J^SJ, is 



A 



(96) 



For the numbers relevant to high- Tc cuprates, the value 
of is very small, i.e. A is large, and, therefore, the 
London limit is well fulfilled (see the estimate in the end 
of subsection IVI CI) . 

The vector potential entering Ea. (|95|l should be inter- 
preted as describing magnetic field averaged over a large 
number of supercells. It is, therefore, possible, that, on 
the scale of I, the true magnetic field fluctuates around 
the exponentially decaying penetration profile character- 
ized by A. 

The relationship between Tc and the zero-temperature 
superfluid density (represented as 1/A^) is plotted in 
Fig. 1101 The points in this plot were obtained by fixing 
the value of the interaction constant g and then vary- 
ing Ea (in Case I) or Sb (in Case II) from zero to very 
large values. The Tc-coordinate of each plot point was 
obtained by taking a pair of values {Sa,g) or {sb,g) and 
then solving numerically Ea. ()46|l or 1)631) for Cases I and 
II, respectively. The 1/A^-coordinate was obtained for 
the same values of {ea,g) or {eb,g) by calculating 5"^ ac- 
cording to Eq.ljnSJ or ijMI). The theoretical plot of Fig.lTUI 
is compared with experiments in Section IVIII Dl 
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FIG. 10: (color online) Critical temperature vs. superfluid 
density (represented as l/A'^) at T = 0: solid line — Case 
I; dashed line — Case II. The plot points are calculated as 
described in the text. 



G. Critical case: Sa ~ st ~ 

In the critical case, all the formulas obtained earher 
have very simple form. In particular, the SC transition 
temperature is given by Ea. (|48|l . which is rewritten here 
as 



9 = 



(97) 



Subst ituting Eq.^ and £a = into Eas. H38l39l40l 
I5UI93|I and using Ea. Ht)2l) when necessary, the zero- 
temperature values of several key quantities can be ex- 
pressed as follows: 



£ao = 4Tc; 



£B — CaoTc', 



Ecs = -2CaoTcN; 



and, finally. 



A 



zo 



he 



2eV -KCaoTc 



(99) 



(100) 



(101) 



The tunneling spectrum for the critical case is shown in 
Fig.Hb). 

I now assume that one model layer represents one 
copper-oxide layer of a real compound, which allows me 
to express the condensation energy per one in-plane cop- 
per atom as 



(102) 



Substituting = 90 K, Z = 4aov/2 « 23A and 
Zo = 6 A into Eas. H98l99ll01ll02l) and recalling that 
Cao = 0.958..., I obtain: e^o = 31 mcV, Eb = 7.4meV, 
Uo = 5.4 K and A = 417 nm. 



VII. DISCUSSION 
A. Realistic features of the model 

In this subsection I list the features, which I expect 
will survive the adaptation of the above model to the 
the properties of real materials, if the 2D diagonal stripe 
hypothesis turns out to be correct. 

The basic underlying feature of this model — the po- 
tential background, which localizes both a- and b-states, 
— should survive the generalizations. In particular, the 
actual shape of the AF domains can be quite distorted 
and thus noticeably different from the perfect diamond 
shape drawn in Fig. ^ It is only important that the 
stripes divide the plane into finite AF domains with the 
alternating sign of the AF order parameter. In the model, 
the disorder in the shape of AF domains can be accom- 
modated through the disorder in the values of Ea , £& and 
9- 

Other supposedly realistic features include the coher- 
ence length of the order of I and the decrease of the crit- 
ical temperature with the increase of |ea — £&!, which is 
associated with the pseudogap. 

Finally, the very unconventional translational symme- 
try of the SC order parameter described in Section IVl El 
should also survive generalizations. 



B. Phase diagram 

In this part, I give the description of the phase dia- 
gram of high-Tc cuprates, which is based, in part, on the 
SC model of Section IVII and, in another part, on a few 
facts borrowed from the next Section, where the model 
predictions are compared with experiments. The phase 
diagram itself is shown in Fig. 1111 

I start at doping concentration Xd ~ 0.06. From the 
viewpoint of the hypothesis adopted in this work, this is 
the lowest concentration at which the 2D stripe pattern 
stabilizes by virtue of some unknown energy balance. At 
this concentration, according to the description given in 
Sections HlllVl and lVIl (i)the length of the stripe supercell 
can be estimated as I — l/{2xd) ~ 8 lattice sites; (ii) 
£a > £b (iii) the value of the pseudogap Sa—eb is maximal; 
and (iv) the superconducting transition temperature is 
minimal. 

At higher doping concentrations, the size of the AF 
domains decreases, eventually saturating at I « 4ao\/2. 
Simultaneously, the pseudogap Eq— also decreases. The 
rationale for the latter assertion is that, at some thresh- 
old doping concentration, stripes should disappear, which 
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FIG. 11: Cuprate phase diagram proposed in the text. 



implies that holes should be expelled from them. Accord- 
ing to Ea. (|^ or JHSI), the decrease of the pseudogap is 
accompanied by the increase of the critical temperature. 

The interaction strength g should also change with the 
increase of doping. It is not obvious a priori whether it 
increases or decreases, but, at least, it is likely that the 
value of g never approaches zero, while the value of ea — Sb 
either reaches zero at some critical doping or becomes 
very small. This suggests that the relative change of 
Ea — eh across the SC doping range has stronger effect on 
the observable quantities than the relative change of g. 

If the value of g were independent of the doping con- 
centration, then the maximum of Tc would be achieved 
at a critical doping corresponding to £a — £6 = 0. How- 
ever, because of the presumed dependence of g on the 
doping level, the optimal doping concentration, a^^o, cor- 
responding to the actual maximum of Tc can be slightly 
shifted with respect to the critical one. The compar- 
ison with experiments in Section IVlll Dl indicates that 
the above notion of the critical doping coincides with the 
"critical doping" Xdc ~ 0.19 identified experimentally in 
j^gfgj5ij52^^ As expected, Xdo ~ 0.16 is not much differ- 
ent from Xdc- The fact that Xdo < Xdc suggests that g 
decreases as the doping level increases. 

Below the critical doping, the inequality Ea > Eb 
implies that the real materials can be describable by 
the model either as Case lA or as Case IIB. (See Sec- 
tion for the discussion of this issue). 

Above the critical doping, it is, in principle, possi- 
ble that either Ea — Eb stays equal to zero, while only g 
changes, or Ea — Eb becomes negative, which means that 
the absolute value of the "inverted" pseudogap starts 
growing again and thus additionally suppresses the crit- 
ical temperature. The analysis of experiments in Sec- 
tion ^^^^ favors the inverted pseudogap scenario and, 
given the model choice between Cases IB and IIA, clearly 
points to Case IB. As discussed in Section |3 the stripe 
superstructure may, for a while, remain stable even after 
Eb becomes greater than Ea- It is further possible that 
the SC transition can contribute to the stabilization of 



stripes by lowering the total energy of the stripe phase. 

In the context of the present proposal, stripes should 
exist in the SC phases of both underdoped and overdoped 
cuprates. It is, however, unclear, whether, in overdoped 
cuprates, stripes can be stabilized without the SC tran- 
sition. Negative answer to this question would imply 
that above the SC transition, overdoped cuprates enter 
a stripeless normal state. Otherwise, the normal state of 
overdoped cuprates may still exhibit some kind of stripe 
order. In turn, if overdoped cuprates enter stripeless 
phase simultaneously with the SC transition, then this 
important aspect cannot be captured by the model of 
Section IVII In such a case, the model scenario for the 
overdoped situations becomes doubtful. 

The conceptual difference between the above descrip- 
tion and the popular idea of competing orders is that, in 
the present proposal, the two orders are not competing, 
but, on the contrary, cooperating: the stripe order is cru- 
cial for the existence of the SC transition, while the SC 
transition can also help stabilizing the stripe order. 



VIII. COMPARISON WITH EXPERIMENTS 
A. Qualitative aspects 

The primary concern in the context of the present pro- 
posal is that no evidence of stripes has been observed so 
far in most optimally doped and overdoped materials. It 
should be pointed out, however, that this proposal stipu- 
lates that strong transverse fluctuations of stripes medi- 
ate superconductivity. If true, this would imply that any 
attempt to observe stripes by pinning charges or freez- 
ing spins, would suppress their transverse fluctuations, 
and thus suppress the SC transition or, at least, signif- 
icantly reduce the critical temperature. Such an inverse 
relation between the amplitude of the stripe fluctuations 
and Tc can explain why the stripes are best observable in 
cuprates having lower such as underdope d^^i*^^ , Nd- 
doped^ and Zn-doped^S-^S*^ LSCO. The transverse fluc- 
tuations of stripes should be strongly coupled to the lat- 
tice. Therefore, some kind of isotope effect should also 
be present in such a system. 

The proper treatment of the single particle excitations 
in the non-superconducting (normal) 2D stripe phase is 
not developed in this work. It is, however, difficult not 
to see that the description of the normal state pseudogap 
given in Section while following only from the basic 
facts about the 2D stripe geometry, bears a strong re- 
semblance of the experimental factsi-^. In particular, the 
disappearance of the pseudogap in the diagonal crystal 
directions can be naturally explained by the presence of 
holes inside the diagonal stripes. In the model frame- 
work, the absolute value of the pseudogap \Ea — Eb\ con- 
stitutes the primary factor suppressing the SC transition 
(see Eas. (|46l63|) '). 

In the SC state, the model predicts such distinctive 
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properties as (i) the suppression of Tc with the growth of 
the pseudogap — ef,| (see Eas. (|46l63ll '): (ii) the emer- 
gence of quasiparticles having coherent dispersion in k- 
space only a,t T < Tc (see the remark fohowing Ea. (|27|l '): 
(iii) the asymmetry of the tunnehng density of states 
(Fig-CK); (iv) hnear density of states around the chemical 
potential (Fig.[7|3,c); and (v) low superfluid density (Sec- 
tion ElFl). In the following subsections HVIIIBIVIIID|I . 
I show that some of the quantitative model predictions 
made without adjustable parameters also agree with ex- 
periments. 



B. Tunneling characteristics 

In this subsection I compare the theoretical tunneling 
spectra at T = with experimental tunneling spectra at 
T ^ Tc- Therefore, the discussion will imply the zero- 
temperature values of all relevant quantities. 

The model predicts two kinds of contributions to the 
tunneling density of states corresponding to A- and B- 
Bogoliubov quasiparticles (see Fig. [TJ . If the Van Hove 
singularities at ±e^o and the delta- function peaks at ±e b 
exist at all, they should be identifiable in experimental 
data. 

The difficulty now is that the tunneling spectra of high- 
Tc cuprates have, in general, only one prominent feature, 
namely, two SC peaks at the opposite values of the bias. 
The energies of the SC peaks are, usually, denoted as ±A, 
and A is referred to as the SC gap. In the superconduc- 
tivity model of Bardeen, Cooper and Schrieffer (BCS), 
A/Tc = 1.76. In high-Tc cuprates, the reported values 
of A/Tc show significant variations but, typically, fall in 
the range between 1.5 and 7, with 4 being a representa- 
tive value. In the present model, the ratios IsaoI/T^c a-nd 
IsbI/Tc can also vary broadly. Their representative val- 
ues are 4 and 1, respectively (see Section IVI C|) . Having 
observed that A/Tc ~ \£ao\/Tc, I identify the experi- 
mental SC peak with the Van Hove singularity in the 
density of A-states. This identification will later prove 
to be consistent with a number of other qualitative and 
quantitative facts. In the following, the variables le^iol 
and A will refer, respectively, to the theoretical and ex- 
perimental values of the same quantity. 

The identification of B-states is more problematic. The 
evidence for their existence is non-trivial, but largely 
indirect. It is based on the STM observations of 
the checkerboard patters in the local density of states 
(LD OS) of Bi-2212is,i9,20^2i discussed in subsec- 

tion lVIIl"n|l . Another possibly related observation is that 
of asymmetric resonance peaks in YBCO by Derro et 
at The asymmetry and the energy range of those peaks 
agree well with the expectations for the delta-peaks due 
to B-states (see Fig. [7{c)). 

There are two possible explanations why, in general, 
B-states are more difficult to observe experimentally. 

The first explanation is that, in a real system, the on- 
site energies Sb can be distributed. As a result, the spec- 



trum of B-states may become broad and featureless. 

The second explanation is that the matrix elements for 
tunneling into b-states (and hence B-states) can be much 
smaller than those for a-states. This, in turn, can be re- 
lated to the fact that b-states are localized in the narrow 
regions inside the stripes, while a-states spread over the 
AF domains and thus have a broader "interface" with 
the environment. Alternatively, it might happen that b- 
states have exotic quantum numbers, in which case tun- 
neling into them can be suppressed at all. 

In the rest of this subsection, I assume that B-states are 
mostly unobservable, and, unless specified otherwise, the 
tests of the model will amount to the comparison between 
the density of A-states and the experimental tunneling 
spectra. 

I limit the choices to the special Cases I and II defined 
in Section IVl Bl Therefore, the model calculation of the 
density of A-states only requires the knowledge of two 
parameters: g and Sa in Case I, or g and Sb in Case II. 
Using Tc and A as input parameters, I can both discrim- 
inate between Cases I and II, and determine the values 
of Ea and Sb- 

The inequality A /Tc > 4 can appear only in the frame- 
work of Case I, in which case, \ea\ and g should be ob- 
tained numerically from Eas. H4t)l50|) (with ns — 1 and 
eao = A). The opposite inequality, A/Tc < 4, can only 
correspond to Case II, i.e. |eb| and g have to be ob- 
tained from Eas. H60l63l65|l (with ub ~ 0, Ca = Cao and 
£^0 = —A). The situation A/Tc — 4 corresponds to the 
critical case described in Section IVI CI 

After the model Hamiltonian is specified, the follow- 
ing tunneling characteristics can be predicted without 
adjustable parameters: (i) the asymmetry in the density 
of A-states or the absence thereof, (ii) the maximum en- 
ergy of A-states sai, and also (iii) the expected positions 
±eB of the delta-peaks representing B-states. 

The selection of Case I implies a strong qualitative pre- 
diction of asymmetric SC peaks. However, the knowledge 
of A and Tc alone cannot help to discriminate between 
Cases lA and IB, which correspond to the opposite "po- 
larities" of the SC peak asymmetry . However, if Case I is 
identified in underdoped cuprates, then the strong expec- 
tation is that Ea > £b, which implies Case lA (larger SC 
peak at the negative bias). The identification of Case II 
implies that the SC peaks are symmetric. (The inequal- 
ity Ea > Eb would then favor Case IIB for underdoped 
cuprates.) 

In Fig. El the model calculations are compared with 
two particularly well resolved STM spectra of Bi-2212. 
The experimental spectrum in Fig. I12f a) was extracted 
from Fig. 2 of RefiM. It is representative of "regular" 
parts of the sample surface (i.e. vortex free and impurity 
free). In this case, the experimentally determined num- 
bers A = 32 meV, Tc = 87 K and A/Tc = 4.3 imply 
Case I with \Ea\ = 7.3 meV and 17 = 31 meV. The the- 
oretical spectrum presented in Fig. Ela) corresponds to 
Ea > 0, i.e. to Case lA. 
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FIG. 12: (color online) Comparison between zero- 
temperature model calculations and low-temperature STM 
spectra of Bi-2212. Solid lines represent the theoretical den- 
sity of states in the same way as in Fig. Dashed lines 
represent experimental spectra extracted from the following 
references: (a) Fig. 2 of Refi— (regular part of Bi-2212 sur- 
face); (b) Fig. 7 of Refi^. 



mental counterparts in the form of the shoulder-hke fea- 
tures located approximately at energies ±e^i predicted 
theoretically. Finally, the delta-peaks shown in Fig. ll2f a) 
are located at ±£b = ±7.3 meV. These are precisely the 
energies, at which, in Ref.-^^ the vortex cores have shown 
anomalous symmetric "humps" absent in the regular SC 
regions. A subsequent study- ^ has revealed that 

The local density of states (LDOS) associated with the 
above humps was later shown to exhibit a checkerboard 
pattern inside the vortex coresi^. In subsection I VIII Cl 1 
will show that this pattern is precisely what one should 
expect from LDOS associated with B-states. 

Now I discuss to what extent the model-based rule 



" A/Tc > 4 <^ asymmetric SC peaks; 
A/Tc < 4 <^ symmetric SC peaks " 



(103) 



is supported by other STM or point-contact 
superconductor — insulator — normal- metal (S-l-N) 
experiments. The values of A used below are obtained 
as half of the difference between the energies of the SC 
peak maxima. 

Supporting evidence: 

Clearly asymmetric SC peaks corresponding to 
A/Tc > 4 have been reported for bi-layer compounds 

Bi-2212 in Refs,5^i-'^6,57,58,59,60,61,62^ 

HgBaaCaCuaOe+a (Hg-I2I2) in Ref&; 

and tri-layer compounds 

BizSrzCasCuaOio+a (Bi-2223) in Ref.^^, 

HgBasCaaCusOg+i (Hg-1223) in Ref.63.65. 

The values of A/Tc extracted from the above references 

cover the range between 4.1 and 6.9 . 

Symmetric SC peaks corresponding to A/Tc < 4 have 
been reported for single layer compounds 
HgBa2Cu04+5 (Hg-1201) in RefiS^, 
TlaBasCuOe (Tl-2201) in Ref.^; 
and bi-layer compounds 

Bi-2212 in Refs.55-56^ also Pb-doped Bi-2212 in RefSL, 
YBCO in Ref.68, 

TlaBasCaCuzO:, (Tl-2212) in Ref.^^ 

The values of A/Tc extracted from the above references 

cover the range between 1.7 and 3.9. 



The experimental spectrum in Fig. ll2f b) was extracted 
from Fig. 7 of Ref.^^*. The relevant experimental numbers 
A = 28 mcV, Tc = 92.3 K and A/Tc = 3.5 imply Case 
II with \eb\ = 4.3 meV and g = 32 meV. The asymmetry 
of the theoretical delta-peaks in Fig. I12r b) corresponds 
to Eb < 0, i.e. to Case IIB. The vertical scale of the 
theoretical plots in Fig. El was chosen to fit best the 
experimental data. 

It thus appears, that not only the model gives the cor- 
rect prediction of the presence or the absence of the SC 
peak asymmetry, but also it predicts the degree of asym- 
metry quantitatively. One should also note that the ter- 
mination points of the theoretical spectra have experi- 



Contradicting evidence: 

Symmetric SC peaks corresponding to A/Tc > 4 have 
been reported for 
Bi-2212 in Refsi55i59i6Qi70. 

Asymmetric SC peaks corresponding to A/Tc < 4 have 
been reported for 
LSCO in Ref-ii, 
Tl-2201 in Ref.S 
Tl-2212 in Ref^. 

However, in the case of LSCO, the difference between the 
peak heights was certainly within the limits of experimen- 
tal uncertainty. In the case of Tl-2201 and Tl-2212, the 
overall impression from the cited references is that the 
SC peaks are largely symmetric. (Most of the spectra 
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reported in the same references and measured on sim- 
ilar junctions pass as symmetric and contribute to the 
"supporting evidence." ) 

Unhke the two spectra shown in Fig. 1121 most of 
the measured tunnehng spectra have more rounded SC 
peaks, which may be the consequence of hmited experi- 
mental resolution. Since a significant broadening of a SC 
peak also shifts the position of its maximum, it is possi- 
ble that, a measured symmetric spectrum indicates the 
ratio A/Tc greater than 4, while the true ratio A/Tc is 
slightly smaller than 4. (One such an example is given 
in Ref.'^^.) The resolution- limited broadening of the SC 
peaks can, therefore, be responsible for, at least a part 
of the "contradicting evidence" in Bi-2212. 

Taken as a whole, the above review of experimental 
data clearly supports the rule pU3l) . Furthermore, this 
rule H1U3|) seems to unify the experimental data, which, 
otherwise, may appear contradicting to each other. (I 
will return to this issue in SectionEmU) 

Another interesting fact is that, in the references cited 
above, the SC peak asymmetry of the bi-layer compounds 
is opposite to that of the tri-layer compounds: the bi- 
layer compounds have the higher SC peak mostly at the 
negative bias (as in Fig. ^^a)), while the tri-layer com- 
pounds have the higher peak at the positive bias. 

Rule can be compared with a more simple pre- 

diction made by Altman and Auerbach^^, that the asym- 
metry of the kind shown in Fig. I12r a') is inherent in all 
high-Tc cuprates. The lack of the asymmetry in some 
of the tunneling spectra and the opposite asymmetry of 
the tri-layer compounds would contradict to the above 
prediction. 

Finally, one additional clear prediction of the model 
is that the inequality A/T^ > 4 implies that, as T ap- 
proaches Tc, the energy of the SC peaks approaches the 
finite value \ea\ (see Eq.ljSUll). The inequality A/Tc < 4 
implies the zero energy of the SC peaks at T — T(, (see 
Ea. (|65|l '). The above prediction is difficult to test, be- 
cause the SC peaks tend to be totally "washed out" in 
the vicinity of Tc. Nevertheless, one can observe that the 
first part of this prediction is consistent with the trend 
in the tunneling data from Refsi^2i22iSi24iS,. The sec- 
ond part is more difficult to test, but it also appears to 
be consistent with the results reported in Refs.^^i£Si21, 
though the results from Refsi2i2^ leave either ambiguous 
or the opposite impression. 



C. B-states and the checkerboard patterns 
observed by STM 

The only experimental evidence, which, at the mo- 
ment, I can identify with B-states is the checkerboard 
modulation of the LDOS observed by STM in the vor- 
tex cores^®. The analogous modulations observed in the 
normal state of Bi-2212'^^ can be attributed to b-states. 
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FIG. 13: (color online) Checkerboard pattern formed by the 
centers of stripe elements. Circles indicate the regions, where 
the densities of B-states in the SC phase or b-states in the 
normal phase are expected to be particularly high. 



(For the alternative interpretations of these (and other) 
checkerboard patterns see Refsi"ii9,20,3i,78,79,80,8i,82)^ 

The checkerboard modulations revealed by the above 
experiments have periodicity of approximately 4 lattice 
periods along the principal lattice directions. From the 
viewpoint of the present proposal, this periodicity is due 
to the response from B-states in the SC state or b- 
states in the normal state. These model states are lo- 
calized around the centers of the stripe elements, which, 
as shown in Fig. 1131 form a pattern with the required 
orientation of the "checkerboard" . 

As a corroborating evidence, I can mention that, in 
the vortex core experimentsi^i^, the LDOS modulation 
was most pronounced at energies approximately equal 
to ±7 meV, in almost exact correspondence with ±£b 
obtained in Section IVlII Bl iust from the knowledge of A 
and Tc for the spectrum presented in Fig. I12f a). In the 
normal state experimental, the typical range of energies, 
where the modulations were most pronounced, was also 
consistent with the possible energies of b-states. 

In regular SC regions of nearly optimally doped B- 
2212, a different kind of LDOS modulations showing 
strong energy dependence has also been observedi^iSS. 
From the view point of the present proposal, this en- 
ergy dependence can be caused by two factors: (i) the 
crossover from the pattern corresponding to B-states to 
the pattern corresponding to A-states; and (ii) the defect- 
induced interference of A-quasiparticles. The first factor 
can be appreciated after one observes that the spatial 
patterns of A-states and B-states are characterized by 
different sets of wave vectors. The pattern of A-states 
has diagonal periods I coinciding with that of the un- 
derlying stripe superstructure, while the periods of B- 
states (shown in Fig I13|l are equal to I / -\/2 and oriented 
along the principal lattice directions. Therefore, as the 
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energy probed by STM increases, the pattern represen- 
tative of B-states gradually transforms into the pattern 
representative of A-states, hence the energy dependence 
of the characteristic wave vectors. The description of 
the second (interference) factor would require quantita- 
tive analysis extending beyond the scope of the present 
work. (The idea, that a different kind of the quasiparti- 
cle interference can entirely explain the energy-dependent 
modulation patterns, was advocated in Refsi^MS. 

Since, in the 2D diagonal stripe picture, the exper- 
imentally observed checkerboard periodicity of 4ao im- 
plies the true underlying period I — AaoV^ along the 
diagonal directions, a direct test for the existence of the 
diagonal superstructure can consist of reconstructing the 
position of hypothetical diagonal stripes from the knowl- 
edge of the checkerboard pattern at lower energies, and 
then checking whether the LDOS modulation at higher 
energies has more pronounced correlations between the 
(approximately) equivalent positions belonging to differ- 
ent supercells. 

Finally, the in-stripe hole content corresponding to the 
present interpretation can be estimated by substituting 
Xd ~ 0.16 and / = aoV^/l = 1/4, into Eq.JH), which 
gives c = Xd/f = 0.6. This number is only slightly 
greater than 0.5 extracted in Section ^ from the ex- 
perimentally observed INS peak splitting in underdoped 
LSCO compounds. 



D. Superfluid density 

In the model framework, the calculation of the critical 
temperature Tc and the superfluid density, S^, requires 
the knowledge of three numbers: \ea ~ g and the 
prefactor of S^. In general, one should expect that both 
|£a — Ebl and g change as functions of doping concentra- 
tion. However, since \ea — eb\, presumably, approaches 
zero not far from the optimal doping, the relative ef- 
fect of this change on the observable quantities should be 
stronger than the effect of the change of g. Therefore, for 
a given family of high-Tj, cuprates, one can obtain an ap- 
proximate relation between 5*0 and Tc by fixing the value 
of g and then varying in Case I or Sb in Case II. The 
two theoretical curves shown in Fig. 1101 were obtained 
precisely in this way. 

In this subsection, I test the model relationship be- 
tween and Tc by superimposing the (rescaled) the- 
oretical plot of Fig. ^1 on the experimental results 
for Tl-220lSmS^, Tlo.s-yPbo.s+vSrzCai-.Y.CusOr 
(Tl-1212)a, Rg-UOM, LSC053,87^ Bi-2212S, Ca-doped 
YBCO (Y:Ca-123)53 ^nd YBCO^^, which report either 
relaxation rate a measured by muon spin rotation (/iSR) 
technique, or the inverse square of the penetration depth 
A extracted from the /j,SR data, field-dependent thermo- 
dynamic measurements, or electron spin resonance (ESR) 
studies. Both a and A^^ should be proportional to 3^,- 
The result is shown in Figs. 1141 and 1151 

In each of Figs. I14r b.d.e.f'). the theoretical plot is 
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FIG. 14: (color online) Critical temperature vs. superfluid 
density (presented as a or A"^ ) at T ^ 0. The theoretical 
plots (solid and dashed lines) are obtained by simple rescaling 
of the plot presented in Fig. IIUI The experimental points are 
extracted from the following references: (a) filled circles - 
ReL% open circles - Ref.^, squares - RefiS; (b) Refi^; (c) 
Ref,-; (d) fiUed circles - Ref-^, open circles - ReffS; (e,f) 
Ref^; (g,h) Refi^. The doping ranges corresponding to the 
experimental points are shown in Fig. 1151 



24 



UD OD 




FIG. 15: (color online) All the data points from Fig. ^| 
rescaled in such a way that the theoretical plots coincide with 
each other. The arrows indicate the direction, in which the 
doping concentrations increase. "UD" and "OD" indicate un- 
derdoped and overdoped samples, respectively. 

rescaled in such a way that the critical case point (the 
one, where the the dashed and the solid curves intersect) 
coincides with the experimental critical point. The latter 
is defined as the point, where the derivative of the experi- 
mental Tc vs. 5*0 dependence undergoes an apparent dis- 
continuous change. In Refsi^iiSiS, the same point was 
found to correspond to the so-called "critical doping con- 
centration" approximately equal to 0.19. In Figs. lT^ g.h) 
the theoretical critical point simply matched the experi- 
mental point of the maximal superfluid density. (Here I 
ignored the issue of anisotropy and treated a-axis and fa- 
axis data as independent data sets.) In Figs. ll4r a.c). the 
rescaling relied on the overall best fit to the experimental 
data rather then on matching a specific point. 

The comparison with experiments in Fig. 1141 amounts 
to the following (crude) quantitative test of the model: 
After the interaction constant g and the prefactor of 
are fixed by establishing the absolute scale for the theo- 
retical critical point, the values of two independent num- 
bers and Tc for non-critical points are obtained by 
varying only one theoretical parameter: Sa in Case I or 
Eb in Case II. 

As evident from Figs. I14land ll5l the theoretical "fish- 
like" plot clearly captures the main features of the ex- 
perimental data, namely: the existence of two different 
regimes, with a critical point of the maximal superfluid 
density separating them. The quality of the quantita- 
tive agreement in Figs. IT^ a.f) is particularly surprising, 
given the crudeness of the test and the fact that the data 
extend into the overdoped region, where the model as- 
sumptions appear less reliable a priori. 

At the point corresponding to the critical case, one can 



also estimate the absolute value of the penetration depth 
from Ec ljlUlll . This estimate, which only requires the 
knowledge of the critical temperature and the transverse 
distance per one Cu02 plane, was made in Section IVI CI 
with the numbers close to those of YBCO or Bi-2212. 
The number obtained (A = 417 nm) is about 3-4 times 
greater than the numbers typically cited for YBCO (see, 
e.g., ReffSi,) and about 2 times greater than the numbers 
cited for Bi-2212 (see, e.g., RefiS). This comparison is 
representative of the general trend: the theoretical for- 
mula H1(J1|I overestimates the penetration depth by about 
factor of three. 

For a simple estimate, which involves only the funda- 
mental constants and two well-known material parame- 
ters {Tc and zq), the factor-of-three agreement with the 
experimental numbers is quite reasonable. One should 
also be conscious of the possibility (mentioned in Sec- 
tion lVI F|) that the profile of magnetic field may strongly 
fluctuate within the penetration depth layer. At the same 
time, the absolute values of the penetration depth are 
typically extracted from the experimental data on the 
basis of theoretical formulas, which do not take into ac- 
count such a possibility. 



E. Correspondence between the model regimes and 
the doping concentrations 

In this subsection, I attempt to establish the corre- 
spondence between model Cases lA, IB, IIA and IIB and 
the doping concentrations of high-Tc cuprates. On the 
basis of the content of Sections IVl lVlllBI and I V 111 PI the 
following four criteria discriminating between the four 
model Cases can be proposed: 

Criterion 1: The inequality A/Tc > 4 indicates Case I, 
while A/Tc < 4 indicates Case II. 

Criterion 2: When the theoretical "fish" -plot from 
Fig. llUl is superimposed on the experimental dependence 
of Tc on the superfluid density, the proximity of a data 
point to the solid line indicates Case I, while the proxim- 
ity to the dashed line indicates Case II. 

Criterion 3: The asymmetry in the tunneling density 
of states characterized by a larger SC peak at negative 
voltages indicates Case lA. The opposite asymmetry in- 
dicates Case IB. 

Criterion 4- When the asymmetry in the tunneling 
density of states is not accessible, I will rely on the postu- 
late, that, in underdoped cuprates, Sa > Sb, which favors 
Case I A over IB, and IIB over IIA. 

Criterion 2, when applied to Fig. 1151 suggests a very 
simple picture: The cuprates are describable by Case II 
at subcritical doping concentrations Xd < Xdc ~ 0.19, 
and by Case I at the supercritical concentrations Xd > 
Xdc ■ Criterion 4 then further narrows the choice to Case 
IIB for Xd < Xdc- This identification implies that, at low 
doping, ^ — Sa > £b, and then, as the doping concen- 
tration increases, the difference Sa — £b decreases until, 
at the critical doping, it becomes equal to zero. The 
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identification of Case I for Xd > Xdc then suggests that 
ka — £b\ starts increasing again (but now with ^ = e^)- 
Since the derivative of function of doping is 

unUkely to change sign exactly at Xd = XdC: I conclude 
that the model pseudogap Sa — £b changes sign as Xd 
passes Xdc- This means that, at the supercritical doping 
concentrations, Sa < Sb, which implies Case IB. Thus the 
assignment following from the above discussion is: 

Xd < Xdc ^ Case IIB . 
Xd > Xdc =^ Case IB. ^ ' 

The clear systematics of the superfluid density data 
should now be contrasted with a less systematic picture 
emerging from the tunneling data. 

Most of the point-contact and STM tunneling spectra 
discussed in subsection IVIII Bl as well as break junction 
and interlayer tunneling spectra, are collected at dop- 
ing concentrations Xd < Xdc- In this doping range, 
the tunneling data for Tl-2201^, Tl2Ba2CaCu208-H5 
(Tl-2212)69, Hg-1201iaS4, LSCOii, YBCOSiaffi (and 
also YbBa2Cu307-afS2i2a) support inequality A/Tc < 4, 
which, according to Criterion 1, suggests Case II in agree- 
ment with the assignment (|1U4|) . 

The tunneling studies of Bi-2212 do not reveal a coher- 
ent picture either in terms of the ratio A/T^ or in terms 
of the SC peak asymmetry. Most of more recent Bi-2212 
tunneling data for Xd < Xrfc^^'^^i''^^'^°i^''^°i^''^^i^^'^' show 
A/Tc > 4, and whenever the asymmetry is evident in the 
data, it mostly points to Case lA — in clear contradiction 
with the assignment (|104|l . At the same time, many other 
(and some of the same) tunneling studiesS^i^S*S£i24i2Si2^ 
find in the same doping range the gap values correspond- 
ing A/Tc ^ 4 and thus, according to Criterion 1 support 
the assignment (|104ll . 

The remarkable fact is that, at least on two 
occasions^SiSiii^ the tunneling spectra of Bi-2212 pre- 
sented in the same paper and measured on samples with 
nearly equal critical temperatures have shown two differ- 
ent values of A/Tc — one significantly greater than 4, 
and the other one smaller than 4. 

The tunneling phenomenology of Bi-2212 can be ex- 
plained by the existence of two different SC states, which, 
for Xd < Xdc, correspond either to Case lA or to Case 
IIB. As argued in Section FV I Ul both SC states are char- 
acterized by sharp minima of the total energy of the sys- 
tem and by the same critical temperature. One of them 
can, e.g., constitute a stable bulk state, while the other 
one a stable or metastable surface state. 

The asymmetric STM spectrum of IIg-1212 reported in 
RefjiS clearly suggest Case lA, which cannot be placed 
within the assignment H104() . 

Criteria 1 and 5, when applied to the asymmetric 
STM spectra for the tri-layer compounds Bi-22232i and 
Hg-1223^^< suggest Case IB, and thus could be compat- 
ible to the assignment 1] 104(1 provided that the doping 
concentration in those samples exceeds the (unknown) 
critical concentration for the corresponding families of 
cuprates. However, the recent interlayer tunneling results 



for Bi-2222^- indicate that, as the doping concentration 
increases, the ratio A/Tc decreases from values larger 
than 4 to values smaller than 4, which, in combination 
with the STM data^'^ rather suggests the assignment op- 
posite to (|104|l . namely: Case IB for Xd < Xdc and Case 
IIB for Xd > Xdc- 

Summary of the findings of this subsection: Assign- 
ment (|104|) is consistent with the superfluid density data 
and/or the tunnehng data reported for Tl-2201, Tl-2212, 
Tl-1212, Hg-1201, LSCO and YBCO. For Bi-2212, the 
same assignment is partially supported by experiments, 
while the overall phenomenology rather suggests the oc- 
currence of two different SC states at the same doping 
concentration. The limited tunneling data on Hg-1212, 
Hg-1223 and Bi-2223 appear to contradict to the assign- 
ment (|104|) with varying degrees of certainty. Assign- 
ment (|104f) . if true, implies that the pseudogap changes 
sign at the critical doping concentration. 



F. Symmetry of the SC order parameter 

The prediction, which distinguishes the present pro- 
posal from many others, is the non-trivial translational 
symmetry of the SC order parameter. This property 
has not been yet investigated experimentally. One 
straightforward experimental test would be to construct 
a nanoscale probe, which is sensitive to the SC phase dif- 
ference between two points separated by one period of 
the expected stripe superstructure. Slightly prior to the 
appearance of the present work, the position-dependent 
sign change of the SC order parameter has also been pro- 
posed by AshkenaziSS in the context of a different stripe- 
related model. 

An important question is whether the SC order param- 
eter introduced in this work is consistent with the phase 
sensitive experiments, which are usually interpreted as 
the evidence for d^2_y2 symmetry of a spatially uniform 
SC state (see e.g. RefiSi). Here I am primarily con- 
cerned with the corner SQUID experiments^^ and with 
the tricrystal junction experimentsSSiSi. 

A definitive discussion of these experiments cannot be 
given at this time because of the following two uncertain- 
ties: 

(i) The continuous family of possible SC solutions ob- 
tained in Section rVl CI has not been narrowed to a single 
one. 

(ii) The boundary conditions necessary for the discus- 
sion of the phase sensitive experiments have not been 
specified. 

In view of the above uncertainties, I limit further dis- 
cussion to presenting just one of several possible inter- 
pretations of the well-known tt phase shift observed in 
the corner SQUID experiments^^. This interpretation is 
illustrated in Fig. E| It is based on the same SC phase 
pattern as the one shown in Fig. El^a). The examina- 
tion of Fig. 1161 reveals that the phases "0" and "7r/2", 
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FIG. 16: (color online) Possible geometry of a corner SQUID 
experiment. 



which characterize both and form the same pat- 
tern along both of the SQUID interfaces, which means 
that this aspect of the phase symmetry is unhkely to 
contribute to the relative phase shift between the two in- 
terfaces. However, in addition, the ^pQ-component also 
has the position dependent sign factor, which is indi- 
cated in the centers of the supercells as or "— ". 
The important fact is that, despite the sign change of 
^'a, the Josephson coupling between and the order 
parameter of the conventional superconductor is not av- 
eraged to zero along each of the two interfaces shown in 
Fig. El Now I assume that, at these interfaces, the order 
parameter of the conventional superconductor preferen- 
tially couples to '^a (as opposed to In such a case, 
the opposite signs of ^ a along the two interfaces imply 
the required phase shift of tt. 

The interpretation of the tricrystal experimentsSLSl, 
in which vortices carrying half of the flux quantum were 
observed, is not practical at this time, because it should 
depend critically on the unknown boundary conditions. 
Here, I can only mention, that, the design geometry of the 
tricrystal experiments is such that, if the postulated 2D 
stripe superstructure exists in each of the three adjacent 
crystals, then the three interfaces between the intersect- 
ing stripe superstructures are nearly equivalent to each 
other. Therefore, the explanation of these experiments 
will likely amount to showing that each of the junctions 
generates a tt phase shift. 

Finally, I discuss the observation of the c-axis Joseph- 
son pair tunneling between YBCO single crystals and the 
films of Pb2&. (Pb is a conventional s-wave superconduc- 
tor.) This experimental fact is difficult to explain on the 
basis of the "d-wave" pictureSi. At the same time, in the 
framework of the present proposal, it can be easily inter- 
preted as follows: When a conventional superconductor 



is placed on the top of the 2D striped system shown in 
Fig. El the contribution of the 5'a-component to the 
Josephson coupling changes sign and thus averages to 
zero, but the contribution from the 5't,-component, which 
does not change sign, can lead to a finite Josephson cur- 
rent. 



IX. CONCLUSIONS 

In the present work, I have analysed several conse- 
quences of the hypothesis that holes doped into high- Tc 
cuprates organize themselves in two-dimensional arrays 
of deep stripes. In particular, on the basis of this hypoth- 
esis, I have formulated and solved a model of supercon- 
ductivity. From that model, I have obtained the tunnel- 
ing spectrum and the superfluid density, which show good 
agreement with experiments. The symmetry of the SC 
order parameter derived from the model is different from 
that of d^2_y2 BCS order parameter. The order param- 
eter obtained in this work has two components, at least 
one of which changes sign as a function of the absolute 
position of the pair on the spatial scale of the stripe su- 
perstructure. A number of other features of this proposal 
such as the geometry of the pseudogap and the effect of 
the pseudogap on the superconducting transition tem- 
perature appear to be in qualitative agreement with the 
phenomenology of high-Tc cuprates. The checkerboard 
pattern of LDOS observed by STM has been interpreted 
as coming from the centers of stripe elements in the 2D 
arrangement of diagonal stripes. This work also indicates 
the possibility, that, in underdoped cuprates, there may 
exist two different kinds of SC states, and that, at the 
critical doping concentration, the pseudogap may change 
sign. 

Even if a future work demonstrates the inadequacy of 
the theoretical assumptions of the present one, the sys- 
tematics of the asymmetry in the tunneling spectra dis- 
cussed in Section lVlII Bl and the scaling of the superfluid 
density data shown in Fig. El should retain the status of 
useful empirical facts. 
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APPENDIX A: APPROXIMATE SOLUTION VIA 
A NON-CANONICAL TRANSFORMATION 



In this Appendix, I present an approximate scheme of 
finding the variational ground state of Hamiltonian Q 
using the following Bogoliubov-like non-canonical trans- 
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formation of a-states in real space 



system when 



Qi = uAi + vrj^y^^Aj; 



(Al) 



together with a regular BogoUubov transformation of fa- 
states 



= sB,, 



jB2 



(A2) 
(A3) 



where and Bij ,j are the annihilation operators of the 
Bogoliubov quasiparticles; u, v, s and w are the trans- 
formation coefficients. These coefficients can be chosen 
real. They must then obey the following normalization 
constraints: 



iv^ = 1, 



w'^ = 1. 



(A4) 



(A5) 



(Variables u and v of this Appendix should not be 
confused with functions u(k) and w(k) defined by 
Eqs.l2M- 

The second term on the right-hand-side of Ea. ljAl|) 
changes sign from supercell to supercell following the sign 
of rji. This sign change is necessary to ensure that the 
canonical fcrmionic anticommutation relations between 
operators Ai and Aj corresponding to neighboring AF 
domains are not violated in the first order of v. This 
transformation is still non-canonical, because it violates 
the anticommutation relations in the second and higher 
orders of v. In order to see this, one can assume that 
A-operators represent true fermions and then check the 
anticommutation relation between operators ai and 
corresponding to a pair of next nearest neighbors. 

Despite the fact that transformation (|A1|) is not canon- 
ical, I will substitute it (together with transformations 
(|A2IA3p into the Hamiltonian (0 and then handle A- 
operators as if they were true fcrmionic operators. 

The justification for such a scheme is three-fold: (i) 
The non-canonical transformation (|A1|) is very natural 
for the structure of Hamiltonian 0. (ii) A priori, this 
scheme represents a controllable approximation in the 
case of small v (large (iii) For arbitrary values of 

parameters Sa, et and g, the ground state energy and the 
quasiparticle excitation energies obtained in the present 
framework turn out to be very close to those obtained 
with the help of the fully canonical transformation of 
Section Ivrni 

Transformation (|A1IA3|) minimizes the energy of the 



w 




where 



X = eail~2nA); 

Y = 4eb(l-2nB); 

Z = gil~2nA)il-2nB). 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 
(All) 
(A12) 



Here tia and ub are the occupation numbers of the quasi- 
particle states described by operators A^ and Bij^cr re- 
spectively, i.e. 



1 



UA 



riB = 



exp(^) + 1 
1 

exp(^) + 1' 



(A13) 
(A14) 



In Eas. l|A13IA14|l . variables ea and eb are the energies 
of the respective quasiparticle states. 

At T = 0, the constraints |u| < 1 and \s\ < 1 impose 
the following condition for the existence of the physical 
solution: 



(A15) 



This condition is satisfied for any non-zero value of g, 
when either Ea = oi Et = 0. 

In the rest of this Appendix, I limit the calculations 
only to Cases lA and HA (in the classification of Sec- 
tion IVTHl. 

Case lA: Eb = 0, Ea > 0. 

Condition Et = implies that, according to 
Eas.f HA8IA9p l. s = 1/v^ and w = -1/^2, both inde- 
pendent of temperature. The coefficients u and v given 
by Eas. ljA6IA7|l with Y ^ have temperature-dependent 
values. 

A natural sign convention for this case is > and 
Eb < 0, i.e., at r = 0, n^i = and hb — 1. 
One can then obtain the energies 



and 



EA^Vel+g^{l-2nBy 



g^{l-2nA){l-2nB) 

4:£A 



(A16) 



(A17) 
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of A- and B-quasiparticles, respectively. 

The zero temperature value of ea following from 
Eq. (jA16|l coincides with that of the Van Hove singularity 
(150(1 characterizing the spectrum of the "canonical solu- 
tion". The values of given by Eas. (|A17|) and (|39|l for 
the two solutions are also close to each other. 

The ground state energy in the present case can be 
evaluated as: 



Egs 



N 



(A18) 



If Eqs given by Ea. l|A18|) is compared with the ground 
state energy of the canonical solution fEa. MUl) ). then the 
difference is never greater than five per cent. 

At finite temperatures, in order to obtain n^i, n^, ea 
and Eb one has to solve the system of equations (IA13I 
UT4l rnel ETtI) numerically It is easy to find, however, 
that the above system of equations always has one triv- 
ial solution: ea — Ea, Eb — ^ with ua and ub given by 
Eas. HA13l I18|l . The condition for the existence of the 
second, non-trivial, solution can be found analytically. 
This condition is: T < Tc, where the critical tempera- 
ture Tc is the solution of Eq. 1(46(1 . Thus the remarkable 
fact is that Tc obtained in the framework of the present 
non-canonical scheme reproduces the "canonical" result 
of Section Ivrni 

Case IIA: Ea = 0, £& > 0. 

According to Eas. lA6IA9(l . the condition Ea ^ im- 
plies that u = l/\/2 and v — l/VS, while s and w have 
temperature-dependent values. 

The sign convention in this case is: ea < 0, Eb > 0, 
i.e., at T = 0, ua = 1 and ub ~ 0. The calculation now 
gives: 



g^{l-2nA)il~2nB) , 

^EB 



el + -g^{l-2nA)'; 



Ec. = -2N 



16^ - 



(A19) 
(A20) 

(A21) 



The critical temperature in this case is again the same 
as obtained from the "canonical" solution, i.e. it is given 
by Eq. pil . 

In summary: The non-canonical variational scheme 
based on Eas. ((AllA2IA3(l predicts the same critical tem- 
perature as the canonical scheme of Sectior fVI CI Fur- 
thermore, the non-canonical scheme predicts the ground 
state energy and the important tunneling characteristics 
within a few per cent from the canonical result. The only 
significant feature of the canonical solution missing in the 
non-canonical one is the absence of the gap in the spec- 
trum of A-quasiparticles in Case II. A related conceptual 



detail is that the non-canonical scheme fails to predict 
the coherent dispersion EaO<) of A-quasiparticles. 

I conclude this Appendix with the following comment: 
In the present variational scheme one can easily find 
that neither the variational energy nor the excitation 
spectrum will change, if the phases of transformation 
((A1IA3(I are modified in the following way: 



(A22) 



(A23) 
(A24) 



where phases (pij can be different for different pairs of 
indices i and j. 

The freedom to vary phases ipij in Eas. l(A22IA24p is, 
at least in part, due to the fact that the parameter space 
of non-canonical transformations is larger than that of 
canonical ones. Therefore, one may try to choose phases 
in Eas. ((A22IA24(l such that transformation ((A22(l be- 
comes canonical. I have found, that, in this way (with 
the selection of phases shown later in Fig. , the next- 
nearest-neighbor anticommutation test described earlier 
in this Appendix can, indeed, be satisfied. However, a 
similar test for the pairs of supercells separated by one 
common neighbor cannot be satisfied independently of 
the choice of phases (py. Transformation ((A22() can- 
not be made rigorously canonical, because it involves 
only the nearest neighbors. It is, however, possible to 
conform with the canonical anticommutation relations, 
if transformation ((A22p is modified to include more re- 
mote neighbors. Such a canonical transformation is much 
easier to describe in fc-space — subject of Section FVI CI 
One can thus conclude that the coherent dispersion of A- 
quasiparticles in fc-space is protected by the Fermi statis- 
tics. 

From a different perspective, one can also observe that, 
in the canonical scheme, the non-intuitive combination of 
phases fa given by Eg . 1(421 minimizes the pairing ampli- 
tude between more remote neighbors, which leads to the 
maximum energy gain from the nearest-neighbors inter- 
acting via the Hamiltonian Q. 



APPENDIX B: MINIMIZATION OF ENERGIES 
(lloll AND deill WITH RESPECT TO PHASES 

The minimization procedure presented in this Ap- 
pendix is equally applicable to the total energy expres- 
sions for both Case I (Eq.lgOI)) and Case II ('Ea. (|^ '). 
Below, in order to be specific, I will focus on the ex- 
pression ((40(1 . This expression can be considered as an 
implicit function of phases entering it through the 
dependence on a single function |y(k)p. 
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From Eg. lf^ . 



|F(k)|^ 



4 |cos^ 
+2 cos 



kRi 4 
kRzH 
kRi 
kR, 



2 
2 
2 



X COS 



(Bl) 



I first note that |y(k)p is the periodic function of k in 
the directions of Ri and R2. Replacing sum in Ea. H4U|) by 
the integral according to the prescription H41(l . I further 
note that the symmetry of the function |y(k)p is such, 
that any shift of the integration region does not change 
the value of the integral. Therefore, the result of the 
integration does not depend on the values of '^^^'^^ and 
^£2i^£±. In the following, in order to be specific, I choose 
^ and = 0. 

The only phase combination to be constrained by the 
minimization of energy H4UI) is v^i+^i-^i-v^ ^ Now I 
switch back to the language of summation and note that 
the summation points in Ea. (|40|l can be divided in pairs 
(k, k') as shown in Fig. (|17|l . Point k is chosen inside the 
white area surrounded by the dashed line, while k' is the 
nearest mirror image of k in the dark area outside of the 
dashed line. For each such a pair. 



while 



where 



|y(k)p = 4[Mk)+p(k) c], 



i/(k')|^ = 4[Mk)-p(k) c], 



/i(k) = cos^ [kRi] + cos^ [kRa] ; 
p(k) = 2cos [kRi] cos [kRa] ; 

C = cos 

Thus the energy H4U|) can be presented as 



(B2) 



(B3) 



(B4) 
(B5) 

(B6) 



k 



-p(k)C] + i^[/i(k)~p(k)c 



(B7) 

where symbol in the sum superscript implies that 
the summation is limited to the area shown in Fig. H17|) 
inside the dashed line. Function F is implicitly defined 
by equation 
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FIG. 17: Two examples of the pairs k-points referred to in the 
text: (ki,k'i) and (k2,k2). Thick solid lines represents the 
boundary of the first Brillouin zone of the stripe superstruc- 
ture. The first k-point of each pair should belong to the light 
region inside the dashed boundary. The second point should 
be obtained from the first one by reflection with respect to 
the nearest dashed line. 



where eyi(k) and ^^(k) = n^[e^(k)] are expressed as 
functions of |V'(k)|2 with the help of Eqs. H29I38|) . Even 
without specifying function F explicitly, one can take the 
derivative of E with respect to C, to find: 

hF ~ 

— = -^p(k){i^'[Mk)+p(k)C] 

-i^'[/i(k)-p(k)C]}, (B9) 

where F' is the first derivative of function F . Each term 
in the sum (|B9|1 is equal to zero, when C = 0, which 
implies an extremum of I have examined a large 
number examples numerically and have found that, in 
all cases considered, the above extremum corresponds to 
the global maximum. This result is also easy to derive 
analytically in the critical case by showing that all pairs 
of terms i^[/i(k) + p(k)C] + F[/i(k) - p(k)C] in Ea.dBTll 
simultaneously reach their maximal values, when C = 0. 
Condition C = substituted into Ea. ljB6|l . then gives: 



^A- ^Pl- 



0, 



(BIO) 



from which Ea. (|42l) follows. 



= [1 - 27iA(k)] £A(k) - £a, 



(B8) 



APPENDIX C: DERIVATION OF THE 
EQUATION FOR THE CRITICAL 
TEMPERATURE 

In order to obtain the critical temperature in Case I of 
Section lVl CI I substitute into Ea. H39|) the limiting values 
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of all quantities as T — Tc —> 0_. In this limit, 

£A(k)^£a, (CI) 

n^(k) ^ (C2) 

exp (tT j + 1 

e_B 0, and tib 1/2, i.e. 2nB — 1^0. In order to 
resolve the uncertainty associated with substituting the 
limiting values of and 2nB — 1, it is necessary to keep 
the next order of Eb in the expression for 2nB — 1, i.e. 

2ns -1--:^. (C3) 



The substitution of Eqs.lEU EH ESj into Eq.^ leads 
to the following equation: 

exp ( ^ ) - 1 1 

= — r \ \ 1 M E I^Wp. (C4) 

Wsa [exp (ft) + ij k 

Now I note, that, independently of the choice of phases 
in Ea.l|BT|l. 

^|y(k)p = 27V. (C5) 

k 

The substitution of Ea. (jC5|) into Ea. (|C4|l then gives 
Eq.gH). 
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